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Abstract 

We explore the prospects to control by use of time-dependent fields quantum trans- 
port phenomena in nanoscale systems. In particular, we study for driven conductors 
the electron current and its noise properties. We review recent corresponding theo- 
retical descriptions which are based on Floquet theory. Alternative approaches, as 
well as various limiting approximation schemes are investigated and compared. The 
general theory is subsequently applied to different representative nanoscale devices, 
like the non-adiabatic pumps, molecular gates, molecular quantum ratchets, and 
molecular transistors. Potential applications range from molecular wires under the 
influence of strong laser fields to microwave- irradiated quantum dots. 
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Notation 



n wire site index, n = 1, . . . ,N 

£ = L, R, lead index 

ne wire site attached to lead £: ni, = 1, n-R — N 

a, P Floquet state indices 

k side-band/ Fourier index 

ea + ih'ja complex quasienergy 

Q, driving (angular) frequency 

T = 27r/0, driving period 

ksT Boltzmann constant times temperature 

r^(e) spectral density of lead £ 

E self energy 

\n) wire site, n = 1, . . . ,N 

l^ait)) — Sfe 6xp(— iA;Oi)|Ma,A;), Floquct state for finite self energy 

\ua,k) kth Fourier coefficient of Floquet state \ua{t)) 

Floquet state for self energy E = 

Ea, \a) eigenenergy and cigenstate of a static Hamiltonian 

Pa/3{t) = (c^Cq), single particle density matrix in Floquet basis 

f{x) — [exp{x/kBT) + Fermi function 
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1 Introduction 



As anticipated by Richard Feynman in his visionary lecture "There's plenty 
of room at the bottom" [1], we witness an ongoing progress in the study of 
physical phenomena on ever smaller scales. Partly, this has been made pos- 
sible by the continuous technical achievements in fabrication and miniatur- 
ization of electronic devices. However, it was the invention of scanning probe 
microscopes [2], which brought about the realization of Feynman's dream, 
namely the selective manipulation of matter on the nanoscale. Since then, 
much progress has been made in nano sciences. In particular, the field of molec- 
ular electronics has emerged, which deals with the realization of electronic 
devices based on the properties of a single or a few molecules. The theoretical 
proposal of a molecular rectifier by Aviram and Ratner [3] has been trend- 
setting for investigating the distinct features of electrical transport on the 
nanoscale. On the experimental side, an ancestor of molecular electronics was 
the pioneering work by Mann and Kuhn [4] on transport through hybrid acid- 
salt surface adlayers. The ongoing advance in contacting single molecules by 
nano-electrodes allows one to perform transport measurements [5-9] . In these 
experiments, the quantum nature of the electrons and the quantum coherence 
across the wire, which is connected to adjacent macroscopic lead electrodes, 
influence various physical properties such as the conductance and the corre- 
sponding current noise statistics. The rapid evolution of molecular conduction 
is documented by recent monographs and article collections [10-13]. 

For the corresponding theoretical investigations, two lines of research are 
presently pursued. A first one starts out from the ab-initio computation of 
the orbitals relevant for the motion of excess charges through the molecular 
wire [14-18] At present, however, the results of such computations generally 
differ by more than one order of magnitude from experimental data, possibly 
due to the equilibrium treatment of exchange correlations [19]. The second line 
employs corresponding phenomenological models in order to gain a qualitative 
understanding of the transport mechanisms involved [20-25]. Two particular 
problems addressed within model calculations are the conduction mechanism 
in the presence of electron-phonon couphng [21-23,26-33] and the length de- 
pendence of the current-voltage characteristics [20,24]. The present work also 
employs rather universal models: We describe the molecules by a linear ar- 
rangement of tight-binding levels with the terminating sites attached to leads. 
Still it is possible to suitably parametrize such tight-binding models in order 
to obtain qualitative results for real systems [34-36] . Furthermore, these mod- 
els also capture the physics of the so-called artificial molecules, i.e. coupled 
quantum dots and quantum dot arrays [37, 38] . 

One particular question that arises in this context is the influence of excitations 
by electromagnetic fields and gate voltages on the electron transport. Such 
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excitations bear intriguing phenomena like photon-assisted tunnehng [38-41] 
and the adiabatic [42-44] and non-adiabatic pumping [45, 46] of electrons. 
From a fundamental point of view, these effects are of interest because the 
external fields enable selective electron excitations and allow one to study their 
interplay with the underlying transport mechanism. In practical applications, 
timc-dcpcndcnt effects can be used to control and steer currents in coherent 
conductors. However, such control schemes can be valuable only if they operate 
at tolerable noise levels. Thus, the corresponding current noise is also of equal 
interest. 

An intuitive description of the coherent electron transport through time- 
independent mesoscopic systems is provided by the Landauer scattering for- 
mula [47] and its various generalizations. Both the average current [48-51] 
and the transport noise characteristics [52] can be expressed in terms of the 
quantum transmission coefficients for the respective scattering channels. By 
contrast, the theory for driven quantum transport is less developed. Scattering 
of a single particle by arbitrary time-dependent potentials has been consid- 
ered [53-55] without relating the resulting transmission probabilities to a cur- 
rent between electron reservoirs. Such a relation is indeed non-trivial since the 
driving opens inelastic transport channels and, therefore, in contrast to the 
static case, an ad hoc inclusion of the Pauli principle is no longer unique. This 
gave rise to a discussion about "Pauli blocking factors" [56,57]. In order to 
resolve such conflicts, one should start out from a many-particle description. 
In this spirit, within a Green function approach, a formal solution for the cur- 
rent through a time-dependent conductor has been presented [58,59] without 
taking advantage of the full Floquet theory for the wire and without obtaining 
a "scattering form" for the current in the general driven case. The spectral 
density of the current fluctuations has been derived for the low-frequency ac 
conductance [60,61] and the scattering by a slowly time- dependent poten- 
tial [62]. For arbitrary driving frequencies, the noise can be characterized by 
its zero-frequency component. A remarkable feature of the current noise in 
the presence of time-dependent fields is its dependence on the phase of the 
transmission amplitudes [62-64]. By clear contrast, both the noise in the static 
case [52] and the current in the driven case [63] depend solely on transmission 
prohahilities. 

In Section 3, we derive within a Floquet approach explicit expressions for both 
the current and the noise properties of the electron transport through a driven 
nanoscale conductor under the infiuence of time-dependent forces [63,64]. This 
approach is applicable to arbitrary periodically driven tight-binding systems 
and, in particular, is valid for arbitrary driving strength and extends beyond 
the adiabatic regime. The dynamics of the electrons is solved by integrating the 
Heisenberg equations of motion for the electron creation and annihilation oper- 
ators in terms of the single-particle propagator. For this propagator, in turn, 
we provide a solution within a generalized Floquet approach. Such a treat- 
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ment is valid only for effectively non-interacting electrons, i.e., in the absence 
of strong correlations. Moreover, this Floquet scattering approach cannot be 
generalized straightforwardly to the case with additional electron- vibrational 
coupling. Better suited for this situation is a quantum kinetic equation for- 
malism which, however, is perturbative in both the wire-lead coupling and the 
electron- vibrational coupling [65,66]. 



An experimental starting point for the investigation of the influence of elec- 
tromagnetic fields on molecular conduction is the excitation of electrons to 
higher orbitals of the contacted molecule. In molecular physics, specific exci- 
tations are usually performed with laser fields. The resulting changes of the 
current through a contacted molecule due to the influence of a laser fleld are 
studied in Section 5. In particular, we focus on the modification of the length 
dependence of the conductivity [67,68]. 



An intriguing phenomenon in strongly driven systems is the so-termed ratchet 
or Brownian motor effect [69-74], originally discovered for overdamped clas- 
sical Brownian motion in asymmetric non-equilibrium systems. Counter- 
intuitively to the second law of thermodynamics, one then observes a directed 
transport although none of the acting forces possesses any net bias. This effect 
has been established also within the regime of dissipativc, incoherent quan- 
tum Brownian motion [74-76]. A mesoscopic device related to ratchets is an 
electron pump [42-46, 77, 78] which indeed might be regarded as a localized 
ratchet. Such systems have already been realized in the quantum domain, but 
almost exclusively operating in the regime of incoherent tunneling [79-83] . In 
Section 6, we study the possibilities for molecular wires to act as coherent 
quantum ratchets and explore the crossover from electron pumps to quan- 
tum ratchets. This requires to investigate thoroughly such quantum ratchet 
systems in the coherent tunneling regime [65,84]. 



The tunneling dynamics of a particle in a bistable potential can be altered 
significantly by ac fields. In particular, it is possible to bring tunneling to a 
standstill by the purely coherent influence of a time-periodic driving [85,86]. 
This so-called coherent destruction of tunneling has also been found in other 
systems [87-89]. In Section 7, we address the question whether a related effect 
exists also for the electron transport through a driven conductor between 
two leads. Moreover, we study the noise properties of the resulting transport 
process [63,64,90,91]. 
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1.1 Experimental motivation 
1.1.1 Coupled quantum dots 

The experimental achievement of the coherent couphng of quantum dots [37] 
enabled the measurement of intriguing phenomena in mesoscopic transport 
[38]. A remarkable feature of coupled quantum dots — the so-called artificial 
molecules with the single dots representing the atoms — is that the energy levels 
of each "atom" can be controlled by an appropriate gate voltage. In particular, 
the highest occupied levels of neighboring dots can be tuned into resonance. At 
such resonances, the conductance as a function of the gate voltage exhibits a 
peak. This behavior is modified by the influence of microwave radiation: With 
increasing microwave intensity, the resonance peaks become smaller and side- 
peaks emerge. The distance between the central peak and the side-peaks is 
determined by the frequency of the radiation field which provides evidence for 
photon-assisted tunneling [38-41]. Photon-assisted tunneling through quan- 
tum dots is, in comparison to its counterpart in superconductor-insulator- 
superconductor junctions [92], a potentially richer phenomenon. The reason 
for this is that quantum dots form a multi-barrier structure which permits 
real occupation and resonant tunneling. Therefore, a theoretical description 
requires to also take into account the influence of the field on the dynamics of 
the electrons localized in the central region between the barriers. The quantum 
dot setup used for the observation of photon- assisted tunneling can also be em- 
ployed as an implementation [93] of the theoretically suggested non-adiabatic 
pump [59,94,95]. 

Related experiments have been performed also with single quantum dots ex- 
posed to laser pulses which resonantly couple the highest occupied orbital and 
the lowest unoccupied orbital of the quantum dot [96]. Such a pulse can cre- 
ate an electron-hole pair which in turn is transformed by a transport voltage 
into a current pulse. Depending on their duration, pulses may not only ex- 
cite an electron but also coherently de-excite the electron and thereby reduce 
the resulting current [97]. In the ideal case, the electron-hole pair is excited 
with probability unity and finally yields a dc current consisting of exactly one 
electron per pulse. This effect might be employed for the realization of a cur- 
rent standard. At present, however, the deviations from the ideal value of the 
current are still of the order of a few percent. 



1.1.2 Molecular wires 

During the last years, it became possible to chemisorb organic molecules via 
thiol groups to a metallic gold surface. Thereby a stable contact between the 
molecule and the gold is estabhshed. This enables reproducible measurements 
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of the current not only through artificial but also through real molecules. 
Single molecule conductance can be achieved in essentially two ways: One 
possible setup is an open break junction bridged by a molecule [5,7,98]. This 
setup can be kept stable for several hours. Moreover, it provides evidence 
for single molecule conductance because asymmetries in the current-voltage 
characteristics reflect asymmetries of the molecule [7,99]. Alternatively, one 
can use a gold substrate as a contact and grow a self-assembled monolayer of 
molecules on it. The other contact is provided by a gold cluster on top of a 
scanning tunneling microscope tip which contacts one or a few molecules on 
the substrate [6, 100]. Yet another interesting device is based on the setup of a 
single-molecule chemical field effect transistor in which the current through a 
hybrid-molecular diode is controlled by nanometer-sized charge transfer com- 
plex which is covalently linked to a molecule in a scanning tunneling micro- 
scope junction [101]. Therein, the effect is due to an interface dipole which 
shifts the substrate work function. Naturally, the experimental effort with such 
molecular wires is accompanied by vivid theoretical interest [8, 10,24]. 

Typical energy scales of molecules lie in the infrared regime where most of 
today's lasers work. Hence, lasers represent a natural possibility to excite the 
electrons of the molecular wire and, thus, to study the corresponding changes 
of the conduction properties. At present, such experiments are attempted, 
but still no clear-cut effect has been reported. The molecule-lead contacts 
seem stable even against relatively intense laser fields, but a main problem 
is the exclusion of side effects like, e.g. heating of the break junction which 
might distort the molecule-tip setup and, thus, be responsible for the observed 
enhancement of the conductance [102] . 

In a recent experiment, Yasutomi et al. measured the photocurrent induced 
in a self-assembled monolayer of asymmetric molecules [83]. They have found 
that even the current direction depends on the wavelength of the irradiating 
light. Albeit not a single-molecules experiment, this measurement represents 
a first experimental demonstration of a ratchet-like effect in molecular wires. 



2 Basic concepts 



Before going in medias res and addressing specific quantum transport situa- 
tions, we introduce the reader to our archetypal working model and the main 
theoretical methods and tools. 



9 



2. 1 Model for driven conductor coupled to leads 

The entire setup of our nanoscale system is described by the time-dependent 
Hamiltonian 

H{t) — H^ireit) + iJleads + -f^contacts, (l) 

where the different terms correspond to the wire, the leads, and the wire-lead 
couplings, respectively. We focus on the regime of coherent quantum transport 
where the main physics at work occurs on the wire itself. In doing so, we neglect 
other possible influences originating from driving-induced hot electrons in the 
leads, dissipation on the wire and, as well, electron-electron interaction effects. 
Then, the wire Hamiltonian reads in a tight-binding approximation with N 
orbitals \n) 

n,n' 

For a molecular wire, this constitutes the so-called Hiickel description where 
each site corresponds to one atom. The fermion operators c„, cj^ annihilate and 
create, respectively, an electron in the orbital \n). Note that in the absence 
of driving a diagonalization of the system Hamiltonian would yield the sta- 
tionary eigenvalues of the wire levels. The influence of an externally applied 
ac field with frequency = 27r/T results in a periodic time-dependence of 
the wire Hamiltonian: Hnn'it + T) — Hnni{t). In an experiment, the driving is 
switched on at a specific time and, thus, the Hamiltonian is, strictly speaking, 
not time-periodic. This can be modeled by a slowly time-dependent driving 
amplitude that assumes its ultimate value after a transient stage in the "in- 
finite past" . Within this work, however, we focus on the transport properties 
at asymptotically long times where the amplitude has already settled [86, 103] 
and, thus, the driving can be assumed periodic. This provides the basis for 
the applicability of a Floquet transport theory. 

The leads are modeled by ideal electron gases, 

Pleads = I] eg (cI^Cl, + 4^%^) , (3) 

where (cr^) creates an electron in the state |Lg) (|Rg)) in the left (right) 
lead. The tunneling Hamiltonian 

^contax^ts = (^L«4gCl + ^KqCnqCj^) + h.C. (4) 

estabhshes the contact between the sites |1), |A^) and the respective lead, as 
depicted with Fig. 1. This tunnehng coupling is described by the spectral 
density 

r,(e) = 27r^|l^,,|25(e-e,) (5) 

9 
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Fig. 1. Level structure of a nano-conductor with N = 5 orbitals. The end sites are 
coupled to two leads with chemical potentials /xl and //r = Ml + eV. 

of lead i = L,K which becomes a smooth function if the lead modes are dense. 
If the leads are modeled by a tight-binding lattice, the r^(e) assume a semi- 
elliptic shape, the so-called Newns- Anderson density of states [104], which 
is sometimes employed in the context of molecular conduction [20,105,106]. 
Within the present context, however, we are mainly interested in the influence 
of the driving field on the conductor and not in the details of the coupling to 
the leads. Therefore, we later on often choose for r^(e) a rather generic form 
by assuming that in the relevant regime, it is practically energy-independent, 

r,(e) r,. (6) 



To fully specify the dynamics, we choose as an initial condition for the left 
(right) lead a grand-canonical electron ensemble at temperature T and electro- 
chemical potential //l(r)- Thus, the initial density matrix reads 

Po oc e~^^^'"^'^''~'^^^^~'^^'^^^^'^^'^ , (7) 

where — J2q ^IqCgg is the number of electrons in lead £ and denotes 
the Boltzmann constant multiplied by the temperature. An applied voltage 

V maps to a chemical potential difference /Ur — = with — e being the 
electron charge. Then, at initial time to, the only nontrivial expectation values 
of the lead operators read {cl,g,Ciq) = fi{eq)5(,i'5qq' where /^(e) = (1 + exp[(e — 
IJi()/k-QT])~^ denotes the Fermi function. 
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Below, we specify the wire Hamiltonian as a tight-binding model composed of 
N sites as sketched in Fig. 1. Each orbital is coupled to its nearest neighbor by 
a hopping matrix element A, thus, the single-particle wire Hamiltonian reads 



N-l 

WwireW = -A 5^ + \n+l){n\)+J2[E^ + Xna{t)] \n){n\, (8) 

n=l n 

where stands for the on-site energies of the tight-binding levels. Although 
the theoretical approach derived below is valid for an arbitrary periodically 
driven wire Hamiltonian, we always assume that the time dependence results 
from the coupling to an oscillating dipole field that causes the time-dependent 
level shifts x„a (t) , where Xn — {N+1— 2n) / 2 denotes the scaled position of site 
\n). The energy a{t) — a{t + T) is determined by the electrical field strength 
multiplied by the electron charge and the distance between two neighboring 
sites. 

An applied transport voltage V is mapped to a symmetric shift of the leads' 
chemical potentials, /xr = — /xl = eV/2. Moreover, for the evaluation of the dc 
current and the zero- frequency noise, we restrict ourselves to zero temperature. 
The zero-temperature limit is physically well justified for molecular wires at 
room temperature and for quantum dots at helium temperature since in both 
cases, thermal electron excitations do not play a significant role. 

In a realistic wire molecule, the hopping matrix element A is of the order 
0.1 eV. Thus, a typical wire-lead hopping rate V = O.IA yields a current 
eV/h — 2.56 X 10~^ Ampere and Q ^ lOA/h corresponds to a laser frequency in 
the near infrared, i.e., to wavelengths of the order 1 /xm. For a typical distance 
of SA between two neighboring sites, a driving amplitude A = A is equivalent 
to an electrical field strength of 2 x 10^ V/cm. It has to be emphasized that the 
amplitude A is determined by the local electrical field between the contacts. 
The difference to the incident field can be huge: Model calculations demon- 
strated that the presence of metallic tips enhances the local field by several 
orders of magnitude [107,108]. This explains the observation that the Raman 
scattering intensity increases drastically once the molecules are adsorbed to a 
metallic surface [109,110]. Coupled quantum dots typically [37,38,40] have a 
distance of less than 1 /im while the coupling matrix element A is of the or- 
der of 30 fj,eV which corresponds to a wavelength of roughly 1 cm. The dipole 
approximation inherent to the time-dependent part of the Hamiltonian (8) 
neglects the propagation of the electromagnetic field and, thus, is valid only 
for wavelengths that are much larger than the size of the sample [111]. This 
condition is indeed fulfilled for both applications we have in mind. 
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2.2 AC transport voltage 



Within this work, we focus on models presented in the previous subsection, 
i.e., models where the driving enters solely by means of time-dependent matrix 
elements of the wire Hamiltonian while the leads and the wire-lead couplings 
remain time-independent. However, it is worthwhile to demonstrate that a 
setup with an oscillating external voltage can be mapped by a gauge transfor- 
mation to the model introduced above. Consequently, it is possible to apply 
the formalism derived below also to situations with an oscillating voltage. 

We restrict the discussion to a situation where the electron energies of the 
left lead are modified by an external T-periodic voltage Kc(^) with zero time- 
average, thus 

e^^e^-eV^{t). (9) 

The generalization to a situation where also the levels in the right lead are 
T-periodically time-dependent, is straightforward. Since an externally applied 
voltage causes a potential drop along the wire [112-114], we have to assume 
for consistency that for an ac voltage, the wire Hamiltonian also obeys a 
time-dependence. Ignoring such a time-dependent potential profile enables a 
treatment of the transport problem within the approach of Refs. [115,116]. In 
the general case, however, we have to resort to the approach put forward with 
this work. 

We start out by a gauge transformation of the Hamiltonian (1) with the uni- 
tary operator 

U^{t) = exp I - ict>{t) (clci + ^ cI^Cl,) } (10) 

where 

0(t) = -|/dt'Kc(t') (11) 

describes the phase accumulated from the oscillating voltage. The transfor- 
mation (10) has been constructed such that the new Hamiltonian H{t) = 
Ul^H{t)Uac — ihUl^Ug^^ possesses a time-independent tunnel coupling. Since, 
the operator ci transforms as Ci — > Ci exp(— i0(i)), the matrix elements Hnn'{t) 
of the wire Hamiltonian acquire an additional time-dependence, 

Hr^n'it) ^ H^rr'it) = Hr,n'it)e-"l'^''^^'r.',-Sn,) + eV^{t)5^^5^,^. (12) 

The second term in the Hamiltonian (12) stems from —ihUl^U^. Owing to the 
zero time- average of the voltage Vac(t), the phase (f){t) is T-periodic. Therefore, 
the transformed wire Hamiltonian is also T-pcriodic while the contact and the 
lead contributions are time- independent, thus, H{t) is of the same form as the 
original Hamiltonian (1). 
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2.3 Tien- Gordon theory 



In order to explain the steps in the current-voltage characteristics of 
microwave- irradiated superconductor- insulator-superconductor junctions [92], 
Tien and Gordon [117] proposed a heuristical theoretical treatment which is 
of appealing simplicity but nevertheless captures some essential features of 
driven transport. The central idea of this approach is to model the influence 
of the driving fields by a periodic shift of the energies in the, e.g. left lead 
according €i,q{t) — ej^q + Acos{Qt), cf. Eq. (9). Then the corresponding lead 
eigenstates evolve as 



where Jk denotes the kth order Bessel function of the first kind. The interpreta- 
tion of the Fourier decomposition (14) is that each state consists of sidebands 
whose energies are shifted by multiples of HQ. For the evaluation of the dc 
current, this is equivalent to replacing the Fermi function of the left lead by 



and formally treating the system as time-independent [117]. While this ef- 
fective static treatment indeed captures the photon-assisted dc current, it 
naturally fails to describe any time-dependent response. 

For time-dependent wire-lead models where the driving shifts all wire levels si- 
multaneously, it is possible to map the driving field by a gauge transformation 
to oscillating chemical potentials. Then, the average current can be evaluated 
from an effective electron distribution like the one in Eq. (15) [118-120]. How- 
ever, generally the time-dependent field also infiuences the dynamics of the 
electrons on the wire. In particular, this is the case for the dipole driving (8). 
Then, a treatment beyond Tien-Gordon theory becomes necessary. Deriving 
an approach which is valid in the general case is the objective of Section 3. 

2.4 Scattering approach for static conductors 

In the absence of a driving field, the computation of the coherent transport 
through mesoscopic structures has become a standard procedure [48-51]. The 
crucial idea goes back to Landauer who postulated already in 1957 [47] that in 
the absence of both inelastic effects and electron-electron interaction, conduc- 
tion can be described as a coherent scattering process of independent electrons. 



\Lq)t=exp 




k=—oo 



(15) 



k 



14 



Then, an infinitesimal voltage V causes the current / = GV with the (linear) 
conductance 

G = jT, (16) 

of a one- dimensional conductor, where T is the transmission probability of 
an electron at the Fermi surface. Since conductors may have non-vanishing 
reflection probability 1 —T, the transmission probability does not necessarily 
assume an integer value. The prefactor e^/h — (25.8 kil)"^ is the so-called 
conductance quantum. 

Originally [47], the conductance (16) has been proposed with T replaced by 
T/(l — T). In the beginning of the 1980's, there has been a theoretical debate 
[121-123] whether or not, the reflection coefficient 1 — T has to be included. 
The controversy was resolved by considering four-terminal devices where two 
terminals act as voltage probes and are considered as a part of the mesoscopic 
conductor [124,125]. Then, V represents the probed voltage and the factor 
1/(1 — T) indeed is justified. In a two-terminal device, however, V denotes the 
externally applied voltage and the conductance includes a contact resistance 
and is given by Eq. (16). 

With the same ideas, Landauer theory can be generalized to the case of a 
finite voltage for which the current reads 

I = lJdE [ME) - ME)]t(E), (17) 

with T{E) being the electron transmission probability at energy E. The elec- 
tron distribution in the left (right) lead is given by the Fermi function /l(r) 
with the chemical potential Hl{k) whose difference — = eV is deter- 
mined by the applied voltage. The linearization for small voltages yields the 
conductance (16). The current formula (17) and the conductance (16) have 
been derived from Kubo formula [122-124, 126, 127] and by means of non- 
equilibrium Green function methods [126, 128-130] for various microscopic 
models. In doing so, one usually starts by defining a current operator, e.g. as 
the change of the electron charge eiVL in the left lead, i.e. / = ie[H,Ni\/h. 
Finally, one obtains the expected expression for the current together with the 
relation 

T{E) = ir[G\E) Y.r{E) G{E) Y.i^{E) ] (18) 

between the transmission probability T{E) and the Green function of the 
electrons. The trace sums over all single-particle states of the wire and = 
\n^)^{n(\ denotes the imaginary part of the self-energy of the terminating 
wire sites which results from the coupling to the respective leads. 

In order to obtain an expression for the related current noise, one considers 
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the symmetrized correlation function 



S{t,t')^l{[AI{t),AI{t')U) (19) 

of the current fluctuation operator AI{t) = I{t) — {I{t)), where the anticom- 
mutator = AB + BA ensures hermiticity. For a stationary process, 

the correlation function S{t, t') — S{t — t') is a function of only the time dif- 
ference. Then, the noise strength can be characterized by the zero-frequency 
component 

/oo 
dT5(T), (20) 
-oo 

which obeys S > Q according to the Wiener- Khinchine theorem. In terms of 
the transmission function T{E), the noise strength reads [52] 



S^^^j dEi^T{E)[f^{E)[l-f^{E)] + f^{E)[l-f^{E)] 
+ T{E) [l - T{E)\ [/r(E) - h{E)\^Y 



(21) 



A dimensionless measure for the relative noise strength, is the so-called Fano 
factor [131] 

^=wv '''' 

Note that in a two-terminal device, both the absolute value of the average 
current and the noise strength are independent of the contact I. Historically, 
the zero- frequency noise (20) contains a factor 2, i.e., one considers S' — 2S, 
resulting from a different dcflnition of the Fourier transform. Then, the Fano 
factor is deflned as F = S' /2e\I\. The definition (22) is such that a Poisson 
process corresponds to F = 1. 

The generahzation of the noise expression (21) to driven systems must also 
account for absorption and emission. Owing to this energy non-conserving 
processes, the zero-frequency noise is no longer given solely in terms of trans- 
mission probabilities but also depends on the phases of the transmission am- 
plitudes [62-64]; cf. Eq. (50), below. 



2.5 Master equation 



A different strategy for the computation of stationary currents relies on the 
derivation of a master equation for the dynamics of the wire electrons. There, 
the central idea is to consider the contact Hamiltonian (4) as a perturba- 
tion, while the dynamics of the leads and the wire, including the exter- 
nal driving, is treated exactly. From the Liouville-von Neumann equation 
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ihg{t) = [H{t), g(t)] for the total density operator g(t) one obtains by standard 
techniques [132, 133] the approximate equation of motion 

Q{t) = - i[i?wire(t) + -ffleads, g{t)] 

(23) 

dTfi^contacts) [-f^contacts ~ Tjt), ^'(^)]]- 

The tilde denotes operators in the interaction picture with respect to the 
molecule and the lead Hamiltonian without the molecule-lead coupling, 
X{t,t') = UQ{t,t') X Uo{t,t'), where C/q is the propagator without the cou- 
phng. For the evaluation of Eq. (23) it is essential to use an exact expression 
for the zeroth-order time evolution operator Uo{t,t'). The use of any approxi- 
mation bears the danger of generating artifacts, which, for instance, may lead 
to a violation of fundamental equilibrium properties [134, 135]. 

In order to make practical use of equation (23), one has to trace over the lead 
degrees of freedom and thereby obtains a master equation for the reduced 
density operator of the wire electrons. Subsequently, the reduced density op- 
erator is decomposed into the eigenstates of the wire Hamiltonian i^wire — or 
the corresponding Floquet states if the system is driven. As a further simpli- 
fication, one might neglect off-diagonal matrix elements and, thus, obtain a 
master equation of the Pauli type, i.e., a closed equation for the occupation 
probabilities of the eigenstates [95,136,137]. For driven systems close to de- 
generacies of the quasienergies, however, such a Pauli master equation is not 
reliable as has been exemphfied in Ref . [66] . 



3 Floquet approach to the driven transport problem 

In the following, we present the Floquet approach for our working model of 
Section 2.1. This derivation is rigorous and exact: It is equivalent to an exact 
treatment in terms of a Keldysh Green function calculation [58] . However, the 
chosen Floquet derivation is here more direct and technically less cumbersome. 

We start out from the Heisenberg equations of motion for the annihilation 
operators in lead £, i.e., 

Ciq = ~^^1^i<l ~ '^il ^riv (24) 

where denotes the conductor site attached to lead £, i.e., rii, = \ and 
n^ — N. These equations are straightforwardly integrated to read 

ce,it) = 0^6-'^"^'-'°^/" - f~'"dT e-'^^^/%^{t - r). (25) 
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Inserting (25) into the Heisenberg equations for the wire operators yields in 
the asymptotic hmit to ^ —oo 

i 1 r°° 

CnAt)=-^T.H^e,n'{t)cAt)--J^ dTre{T)Cn,{t-T)+Ut), (26) 

Cnit) =- ^J2Hnn'{t)Cn'{t), 71 = 2, . . . , iV - 1, (27) 
n' 

where the lead response function r^(i) results from the Fourier transformation 
of the spectral density (5), 

r.(^) = /^e-^^*/'^r,(6). (28) 

In the wide-band limit (6), one obtains r^(t) = T£6{t) and, thus, the equa- 
tions of motion for the wire operators are memory-free. The influence of the 
operator-valued Gaussian noise 

Ut) = (29) 

is fully specified by the expectation values 

(Ut)) = 0, (30) 
i^l'it') m) = See I ^ e-^^(*-*')/^ r,(e)/,(6) , (31) 

which for the uncorrelated initial state (7) follow from the definition (29). 
It is convenient to define the Fourier representation of the noise operator, 
Cii^) — I dt exp{iet / h)^j({t) whose correlation function 

(d(e)e£'(e')) = 27rr,(e)/,(e) S{e - e') Sw (32) 
is obtained directly from Eq. (31). 

3. 1 Retarded Green function 

The equations of motion (26) and (27) represent a set of hnear inhomogeneous 
equations and, thus, can be solved with the help of the retarded Green function 
G{t, t') = U{t, t')e{t - t') which obeys 

where T{t) = |l)rL(t)(l| + |A^)rR(t)(iV|. At this stage, it is important to note 
that in the asymptotic limit to — > — oo, the l.h.s. of this equation is periodic 
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in As demonstrated in the Appendix, this has the consequence that for the 
propagator of the homogeneous equations obeys U {t, t') — U{t + T,t' + T) 
and, accordingly, the retarded Green function 

G(t, e) = Hdr e''^'^ U(t, t-T)^G(t + T, e) (34) 
a Jo 

is also T-periodic in the time argument. Thus, we can employ the Fourier 
decomposition G{t,e) = X);- e^''^^* G'*^'^)(e), with the coefficients 

G'('=)(e) =^ r dt e'*^"*G(t,e). (35) 
T Jo 

Physically, G*-'^-' (e) describes the propagation of an electron with initial energy 
e under the absorption (emission) of photons for k > {k < 0). In the 
limiting case of a time-independent situation, G{t, e) becomes independent of 
t and, consequently, identical to G^^^ (e) while all sideband contributions with 
/c 7^ vanish. 

From the definition of the Green function, it can be shown that the solution 
of the Heisenberg equation (26), (27) reads 

roo 

Cn{t)^ih^ drGn,nAt,t-r)Ceit-r). (36) 

Inserting for Gn,,n{t,t') — {n\G{t,t')\n£) the Fourier representation (34), one 
obtains the form 

cn(^) = ^ E / dee-'^'/^Gn,nAt, e) U^), (37) 
which proves more convenient. 

Below, we need for the elimination of back-scattering terms the relation 
G^^it, e') - G{t, e) = (ih^ - e' + ejG^t, e')G{t, e) 



+ i r dTe''^/^'GUt,e')r(T)G(t-T,e) (3J 
Jo 

/■oo 

+ i / dTe-''^/''GUt-T,e')rUT)G(t,e). 
Jo 



A proof of this relation starts from the definition of the Green function, Eq. 
(33). By Fourier transformation with respect to t', we obtain 



■4 

dt 



+ e - n{t)^ G{t, e) + i e''^/" r{r) G{t - r, e) = 1 (39) 



which we multiply by G\t,t) from the left. The difference between the re- 
sulting expression and its hermitian adjoint with e and e' interchanged is 
relation (38). 
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3.2 Current through the driven nano-system 

The (net) current flowing across the contact of lead £ into the conductor is de- 
termined by the negative change of the electron number in lead i multiplied by 

the electron charge — e. Thus, the current operator reads = ie[H{t), N^l/h, 
where = J2q ^IgCig denotes the corresponding electron number. By using 
Eqs. (25) and (29), we obtain 

hit) = 1 dr {re(r)c\(t)c,(t - r) + ri(r)cl(t - r^t)} ^^^^ 
-e{c{it)Ut)+^Kt)ci{t)}. 

This operator-valued expression for the time- dependent current is a conve- 
nient starting point for the evaluation of expectation values like dc current, 
ac current, and current noise. 

3.2.1 Average current 

In order to evaluate the current {Ii,{t)), we insert the solution (37) of the 
Heisenberg equation into the current operator (40) and use the expectation 
values (32). The resulting expression 

(^lW) = X?/^' fdT{e'^-/'^GUt,e)r^iT)Gu{t-T,e)T,{e)Me) 

+ e-'^^/^GUt - T, e) r* (r) Gu{t, e)Te{e)fe{e)) 
+ iejde (Gl,{t,e) - Giiit,e))ri^ie)Me) 

(41) 

still contains back-scattering terms Gu and, thus, is not of a "scattering form" . 
Indeed, bringing (41) into a form that resembles the static current formula 
(17) requires some tedious algebra. Such a derivation has been presented for 
the linear conductance of time-independent systems [126] , for finite voltage in 
the static case for tunneling barriers [128] and mesoscopic conductors [129], a 
wire consisting of levels that couple equally to both leads [58], and for weak 
wire- lead coupling [56]. For the general time-dependent case in the absence 
of electron-electron interactions, such an expression has been derived only 
recently [63,64]. 

Inserting the matrix element (1| . . . |1) of equation (38), eliminates the back- 
scattering terms and we obtain for the time-dependent current the expression 

(hit)) = 11 de {T^t, e)/R(e) - T^^it, e)Me)} - ^^q^(t) (42) 
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where 



QLit) ~ J derj^{e)J2\Gni{t,e)fMe) (43) 

denotes the charge oscillating between the left lead and the wire. Obviously, 
since is time-periodic and bounded, its time derivative cannot contribute 
to the average current. The corresponding charge arising from the right lead, 
qB.{t), is a priori unrelated to qhit); the actual charge on the wire reads 
?l(^)+?r(^)- The time-dependent current is determined by the time-dependent 
transmission probability 

rLR(i,e) = 2Re H dr e'^^/Tj^ir) Gl^it,e) Gi^it - r,e)rnie). (44) 
Jo 

The corresponding expression for TYih{t,e) follows from the replacement 
(L, 1) (R, A^). We emphasize that (42) obeys the form of the current formula 
obtained for a static conductor within a scattering formalism. In particular, 
consistent with Refs. [49,56], no "Pauh blocking factors" (1 — /^) appear in 
our derivation. In contrast to a static situation, this is in the present context 
relevant since for a driven system generally 

7kL(t,e)7^TLR(t,e) (45) 

such that a contribution proportional to /L(eg')/R(^?) would not cancel [56,57]. 

In order to obtain an expression for the dc current, we insert for the Green 
function the Fourier representation (35) followed by performing the average 
over time t. Then, the average current becomes 

^~ - X E / {tS(^)/r(^) - 4S(e)/L(6)} , (46) 



where 



fe=— 00 



TSie) =rL(6 + khn)rnie) G\'>ie) \ (47) 



2 



45(6) =rR(6 + khn)T^{e) Gl^i(e) ", (48) 

denote the transmission probabilities for electrons from the right lead, respec- 
tively from the left lead, with initial energy e and final energy e + khQ, i.e., 
the probability for an scattering event under the absorption (emission) oi \k\ 
photons if A; > (A; < 0). 

For a static situation, the transmission probabilities T^^^^e) and T^^{e) are 
identical and contributions with k ^ vanish. Thus, it is possible to write the 
current (46) in the form (17) as a product of a single transmission probability 
T(e), which is independent of the direction, and the difference of the Fermi 
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functions, /R(e) — fhi^)- We emphasize that in the driven case this is no longer 
true. 

3.2.2 Noise power 

Like in the static case, we characterize the noise power by the zero-frequency 
component of the current-current correlation function (19). However, in the 
driven case, Se{t, t') — Se{t + T,t' + T) is still time-dependent. Since it shares 
the time-periodicity of the driving, it is possible to characterize the noise level 
by the zero-frequency component of Si{t,t — r) averaged over the driving 
period, 

} fT f+oo 

Se = 7F dTSe{t,t-T). (49) 

1 Jo J-oo 

It can be shown [64] that for driven two-terminal devices, Se is independent 
of the contact i, i.e., Si^ — S-ei = S. 

We start by writing Si^{t, t — r) with the current operator (40) and insert the 
solution (37) of the Heisenberg equations. We again employ relation (38) and 
finally obtain 

(50) 

We have defined e^'^) — e + khO, and fi = l — fi. It can be shown (cf. Sec- 
tion 3.5.1) that in the undriven limit, the noise power (50) depends solely on 
the transmission probabilities and is given by Eq. (21). In the time-dependent 

case, however, the noise expression (50) cannot be brought into such a con- 
venient form and, thus, generally depends on the phase of the transmission 
amplitude. 

3.2.3 Floquet decomposition 

For energy-independent wire-lead coupling, i.e. in the so-called wide-band limit 

Ti{e) = Ti the lead response function (28) reads Ti{t) = r^5(t). Consequently, 
the integro-differential equation (33) for the Green function becomes a pure 
differential equation. Then, determining the Green function is equivalent to 
computing a complete set of solutions for the equation 
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h 



de{rR(eW)rR(e) 



(fc') 



+ rR(6('=))rL(e) ErL(6feOG 



(fc'-fc)/ 

liV J 



+ same terms with the replacement (L, 1) <-> (R, N). 



where the self-energy 

S= 11)^(11 + 1^)^(^1 (52) 

results from the coupling to the leads. Equation (51) is linear and possesses 
time-dependent, T-periodic coefficients. Thus, following the reasoning of Ap- 
pendix A, it is possible to construct a complete solution with the Floquet 
ansatz 

\iJa{t)) =exp[{-[ea/fi- la)t\\Ua{t)), (53) 
l^aW) =Z!l«a,fc)exp(-iA;fit). (54) 

k 

The so-called Floquet states \ua{t)) obey the time-periodicity of Ttwheit) and 
have been decomposed into a Fourier series. In a Hilbert space that is extended 
by a periodic time coordinate, the so-called Sambe space [138], they obey the 
Floquet eigenvalue equation [139, 140] 

(Hwiieit) - iS - ih^^\Ua(t)) = (Sa - i/l7a) |lia(i)) ■ (55) 

Due to the Brillouin zone structure of the Floquet spectrum [138, 139, 141], it 
is sufficient to compute all eigenvalues of the first Brillouin zone, —Ml/ 2 < 
€a < Ml/2. Since the operator on the l.h.s. of Eq. (55) is non-Hermitian, the 
eigenvalues — ihrfa are generally complex valued and the (right) eigenvectors 
are not mutually orthogonal. Thus, to determine the propagator, we need to 
solve also the adjoint Floquet equation yielding again the same eigenvalues but 
providing the adjoint eigenvectors Iw^it))- It can be shown that the Floquet 
states \ua{t)) together with the adjoint states \u'^{t)) form at equal times a 
complete bi-orthogonal basis: {u^{t)\u0{t)) — 5^^ and Y.a l'^a(0)('^a (01 — 1- 
A proof requires to account for the time-periodicity of the Floquet states since 
the eigenvalue equation (55) holds in a Hilbert space extended by a periodic 
time coordinate [139,142]. For details, see Appendix A. 

For the special case [59] of a wire with N = 2 sites which couple equally 
strong to both leads, i.e., Fl = Fr, the self-energy is proportional to the unity 
matrix. Consequently, the Floquet states \u^{t)) become independent of the 
self-energy. 

Using the Floquet equation (55), it is straightforward to show that with the 
help of the Floquet states \ua{t)) the propagator can be written as 

U{t,t') = ^e-(^«/'^-^«)(*-*') K{t)){u^{t% (56) 

a 

where the sum runs over all Floquet states within one Brillouin zone. Conse- 
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quently, the Fourier coefficients of tlie Green function read 




hJo T Jo 

\Ua,k'+k){uXk' 



(58) 



(57) 



For tfie exact computation of current and noise, we solve numerically the 
Floquet equation (55). With the resulting Floquet states and quasicucrgics, 
we obtain the Green function (35). In the zero temperature limit, the Fermi 
functions in the expressions for the average current (46) and the zero-frequency 
noise (50) become step functions. Then, the remaining energy integrals can be 
performed analytically since the integrands are rational functions. 

3.3 Symmetries 

A system obeys a discrete symmetry if its Hamiltonian is invariant under a 
symmetry operation S = {S~^)~^, i.e, if S~^H{t)S = H{t). Then the corre- 
sponding transition amplitude in position representation fulfills the relation 



such that the corresponding transmission probabilities are identical. The com- 
plex conjugation in Eq. (59) holds if <S includes time inversion [143]; then the 
r.h.s. becomes {x'\U{t' ,t)\x) . If S\x') ^ |a;), relation (59) means that two dif- 
ferent scattering processes occur with the same probability. Correspondingly, 
in a time- dependent transport problem as defined by the Hamiltonian (1), 
the presence of a symmetry implies that two different transport channels have 
equal transmission probability. 

Here, we identify the channel which is related to ^^^(e) given a certain sym- 
metry is present. In particular, we consider systems that are invariant under 
the transformations studied in the Appendix A. 3 which are combinations of 
the transformations 



For the tight-binding model sketched in Fig. 1, the parity operation (60) maps 
the lead states and the wire sites according to 




(59) 



Sp : X ^ —X, 
ST-.t^ -t, 
SG:t^t + T/2. 



(60) 
(61) 
(62) 



«Sp : {Lq, n) ^ {Rq, N + l-n), 



(63) 
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Fig. 2. Transmission of an electron with energy e under the absorption of k photons 
(soUd line) and its symmetry related process (dashed) for (a) time-reversal symme- 
try, (b) time-reversal parity, and (c) generalized parity. The sketched processes occur 
with equal probability. 

where n = 1, . . . ,N labels the wire sites and Lq (Rg) the states in the left 
(right) lead. Both the parity Sp and the time inversion St can be generalized 
by an additional shift of position and time, respectively. Alternatively, one 
can place the origin of the corresponding axis properly. For convenience, we 
choose the latter option. 



It should be mentioned that for the periodic driving considered in this work, 
the system contains a further symmetry, namely the time-translation by a 
full driving period. This has already been taken into account when deriving a 
Floquet transport theory and cannot be exploited further. 



3.3.1 Time-reversal symmetry 

If the Hamiltonian obeys time-reversal symmetry S^, i.e., if H{t) = H{—t), Eq. 
(59) yields {l\U{t,t')\N) = {N\U{-t' , -t)\l). Inserting into the definition of 
the Green function G*^'^)(e), Eqs. (34) and (35), results in the relation G^^^{e) = 
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^Ni (^+ khil), where we have shifted the hmits of the t-integration using the 
relation G{t, e) = G{t + T, e). Thus, the transmission probabihties obey 

Tit^{e)=Ti^'\e + khQ), (64) 
i.e., the scattering processes sketched in Fig. 2a occur with equal probability. 

A time-independent system in the absence of magnetic fields presents a partic- 
ular case of time-reversal symmetry since all transmissions probabilities with 
k ^ vanish and, thus, T^^i^) — ^lr(^) — T{e). 

3.3.2 Time-reversal parity 

Systems driven by a dipole force with purely harmonic time-dependence obey 
the so-called time-reversal parity Stp = S^Sp, i.e., a combination of time- 
reversal symmetry and parity. This of course implies that the static part of 
the Hamiltonian has to obey spatial parity which requires identical wire-lead 
couplings, rL(e) = rR(e). The consequences for the Floquet states are dis- 
cussed in the Appendix A. 3 while here, we derive the consequences for the 
transmission probabilities. 

By the same reasoning as in the case of time-reversal symmetry discussed 
above, but with additionally interchanging left and right, we find G'^^(e) = 
G[^^ (e -|- khfl) which yields equal transmission probabilities for the scattering 
events sketched in Fig. 2b, i.e. 

Ti2{e)=Til'\e + khQ). (65) 

Interestingly, timc-rcvcrsal parity relates two scattering events that both go 
into the same direction. Therefore, relation (65) has no obvious consequence 
for any dc current. Still time-reversal parity entails an intriguing and more 
hidden consequence for non-adiabatic electron pumping by harmonic mixing 
as a function of the wire- lead coupling [65] . We discuss this effect in the context 
of non-adiabatic electron pumping in Section 6.3. 

3. 3. 3 Generalized parity 

A further spatio-temporal symmetry that has an impact on the transmission 
properties is the so-called generalized parity Sqp = S^Sp, i.e., a parity opera- 
tion combined with a time shift by half a driving period. This symmetry also 
explains qualitatively the quasienergy spectra found in the context of driven 
quantum tunneling [85, 86, 144, 145]. 

If the wire-lead Hamiltonian is invariant under iSqp, the time evolution oper- 
ator obeys {l\U{t,t')\N) = (iV| ?7(t + r/2,t' + r/2) 1 1). Inserting into Eq. (34) 
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results in G[j(f{e) — G'Vi(^) ^^d, thus, the scattering events sketched in Fig. 2c 
obey 

4?(e)=TS(6). (66) 

Again, we have shifted the integration hmits by using the time-periodicity of 
the Green function G{t,e). 

3.4 Approximations 

In Section 3.2, expressions for the current and the noise power have been 
derived for a periodic but otherwise arbitrary driving. Within the wide-band 
hmit, both quantities can be expressed in terms of the solutions of the Floquet 
equation (55), i.e., the solution of a non-Hermitian eigenvalue problem in an 
extended Hilbert space. Thus, for large systems, the numerical computation 
of the Floquet states can be rather costly. Moreover, for finite temperatures, 
the energy integration in the expressions (46) and (50) have to be performed 
numerically. Therefore, approximation schemes which allow a more efficient 
computation are of much practical use. 

3.4-1 Weak-coupling limit 

In the limit of a weak wire- lead coupling, i.e., for coupling constants F^ which 
are far lower than all other energy scales of the wire Hamiltonian, it is possible 
to to derive within a master equation approach a closed expression for the dc 
current [65]; cf. Section 4. The corresponding approximation within the present 
Floquet approach is based on treating the self-energy contribution — iS in the 
non-Hermitian Floquet equation (55) as a perturbation. Then, the zeroth order 
of the Floquet equation 

(n^ireit) - ih^) \Mt)) = ellMt)), (67) 

describes the driven wire in the absence of the leads, where {(pait)) = 
J2k^w{~^kflt)\4>a,k) are the "usual" Floquet states with quasienergies e°. In 
the absence of degeneracies the first order correction to the quasienergies is 
—ifijl^ where 

i'a = ll^Y^umut)) (68) 

-^i:\Wa,k)\' + ^Y.\iNM\'- (69) 

Since the first order correction to the Floquet states will contribute to neither 
the current nor the noise, the zeroth-order contribution \ua{t)) = \u^{t)) — 
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already sufficient for the present purpose. Consequently, the trans- 
mission probability (47) assumes the form [64] 

{N\(t>a,k' ) {(t>a,k'+k 1 1) ( 1 1 (t>p,k"+k) {(t>l3,k" I N) 

^ ^ ' - i^l + k'hn + i;i7i)] - i^l + ^"^^^ - i^7^)] ^ ^ 

and T^^{e) accordingly. The transmission probability (70) exhibits for small 
values of Ti sharp peaks at energies e° + k'hfl and + k"hfl with widths /i7^ 
and h'f^. Therefore, the relevant contributions to the sum come from terms for 
which the peaks of both factors coincide and, in the absence of degeneracies 
in the quasienergy spectrum, we keep only terms with 

a = p, k' = k". (71) 

Then provided that 7^^ is small, the fraction in (72) is a Lorentzian and can be 
approximated by 7r5(e — — k'hfl) /h^yl^ yielding the transmission probability 

tS(6) = T^T^Y.T^\Wa,k'+k){MN)\'6{e - e° + k'hn) (72) 

^Til'\e + khfl). (73) 

The last line which follows by substituting k' ^ k' — k, means that the trans- 
mission probabilities in the weak-coupling limit obey the same relation as in 
the case of time-reversal symmetry, cf. Eq. (64) even in the absence of any 
symmetry. 

The energy integration in (46) can now be performed even for finite temper- 
ature and we obtain for the dc current the expression 

/ = f E [Uel + k'hfl) - Mel + khQ)] . (74) 

a,k,k' La + J- Ra 

The coefficients 

TLafc =rL|(l|0a,fc)P, = E ^Lafc , (75) 

k 

rR.ik=rR|(7V|</,„^,)|2, rR,=ErRafe, (76) 

k 

denote the overlap of the A;th sideband \(l)a,k) of the Floquet state \(l)a{t)) with 
the first site and the last site of the wire, respectively. We have used 2%^ = 
^La+^Ra which foUows from (69). Expression (74) can been derived also within 
a rotating- wave approximation of a Floquet master equation approach [65]; 
cf. Sect. 4.3. 

Within the same approximation, we expand the zero-frequency noise (50) to 
lowest-order in Tf After inserting the spectral representation (58) of the Green 
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function, wc again keep only terms with identical Floquet index a and identical 
sideband index k to obtain 

+(rL + rL)rLa./L(e°+A;nQ)} ^^^^ 

+ same terms with the replacement L <-> R. 



Of particular interest for the comparison to the static situation is the limit of 
a large applied voltage such that practically /r = 1 and /l = 0. Then, in Eqs. 
(74) and (77), the sums over the sideband indices k can be carried out such 
that 

^oo^Ef^^^, (78) 

C rLarRa(rLa + T^q) , , 

1^ a U La + J- Raj 

These expressions resemble the corresponding expressions for the transport 
across a static double barrier [52]. If now Fl^ = Fro, for all Floquet states 
\(t)a{i))-, we find F = 1/2. This is in particular the case for systems obey- 
ing reflection symmetry. In the presence of such symmetries, however, the 
existence of exact crossings, i.e. degeneracies, limits the applicability of the 
weak-coupling approximation and a master equation approach (cf. Sect. 4) is 
more appropriate. 



3.^.2 High-frequency limit 

Many effects occurring in driven quantum systems, such as coherent destruc- 
tion of tunneling [85] or current and noise control [63,90], are most pronounced 
for a large excitation frequency Q. Thus, it is particularly interesting to de- 
rive for the present Floquet approach an expansion in terms of l/fi. Thereby, 
the driven system will be approximated by a static system with renormalized 
parameters. Such a perturbation scheme has been developed for two-level sys- 
tems in Ref. [141] and applied to driven tunneling in bistable systems [86] and 
superlattices [145]. For open quantum system, the coupling to the external 
degrees of freedom (e.g., the leads or a heat bath) bears additional complica- 
tions that have been solved heuristically in Ref. [91] by replacing the Fermi 
functions by effective electron distributions. In the following, we present a 
rigorous derivation of this approach based on a perturbation theory for the 
Floquet equation (55). 

We assume a driving that leaves all off-diagonal matrix elements of the 
wire Hamiltonian time-independent while the tight-binding levels undergo a 
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position-dependent, time-periodic driving /n(t) = fn{t + T) with zero time- 
average. Then, the wire Hamiltonian is of the form 

n^Ut) = ^0 + E/n(^) \n){n\. (80) 

n 

If Ml represents the largest energy scale of the problem, we can in the Flo- 
quet equation (55) treat the static part of the Hamiltonian as a perturbation. 
Correspondingly, the eigenfunctions of the operator J2n fn{'t)\'>T') {iT'I ~ ihd/dt 
determine the zeroth order Floquet states 

e-^^"(*) \n). (81) 

We have defined the accumulated phase 

Fn{t) fnit') = F„it + T), (82) 

which is T-periodic due to the zero time-average of fn{t). As a consequence 
of this periodicity, to zeroth order the quasienergies are zero (mod/lQ) and 
the Floquet spectrum is given by multiples of the photon energy, khil. Each 
k — 0,±1,±2, . . . defines a degenerate subspace of the extended Hilbert space. 
If now hQ is larger than all other energy scales, the first order correction to 
the Floquet states and the quasienergies can be calculated by diagonalizing 
the perturbation in the subspace defined hy k — 0. Thus, we have to solve the 
time-independent eigenvalue equation 

(Heff-iE)|a) = (6i-i%^)|a). (83) 

The static effective Hamiltonian Hes is defined by the matrix elements of the 
original static Hamiltonian Ho with the zeroth order Floquet states (81), 

(Kfr)nn' = / ^ ci^"W(7^o)nn' e'^^-'^ . (84) 

Jo 1 

The t-integration constitutes the inner product in the Hilbert space extended 
by a periodic time coordinate [138] (for details, see Appendix A. 2). To first 
order in 1/^2, the quasienergies e\ — ihry^ are given by the eigenvalues of the 
static equation (83) and, consequently, the corresponding Floquet states read 

\u^{t)) = ^e-^^"(*) \n){n\a). (85) 

n 

The fact that all Fn{t) are T-periodic, allows one to write in (85) the time- 
dependent phase factor as a Fourier series, 

e-iF4*)^^a„,,e-^^*. (86) 

k 
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Thus, {n\ua,k) — o,n,k{'n\a) and the Green function for the high-frequency 
driving reads 

Gin'(e) = E«n,.'+.<^.^Gf„,(e - k'hn), (87) 

k' 

where ^^^(e) denotes the Green function corresponding to the static Hamilto- 
nian Hes with the self-energy E. Finally, substituting e — > e -|- k'HQ and using 
the sum rule J2k' '^n,k+k'0'n,k' — ^k,o, we obtain 

/ = I y de reff(e){/R,eff(6) - /L,efr(e)}. (88) 

The effective transmission probability Teff(e) = Ti^TYi\G'i%{e)\'^ is computed 
from the effective Hamiltonian (84); the electron distribution is given by 

/L,eff(e) = E M'M^ + km) (89) 

k 

and /R,eff follows from the replacement (1, L) (A?", R). 

In order to derive a high-frequency approximation for the zero-frequency 

noise S, we insert (87) into (50) and neglect products of the type G'^^{e — 
khVt) G^^{t — k'MT) for k ^ k'. Employing the above sum rule for the Fourier 
coefficients a„^fe, we obtain for the noise the static expression (21), but with 
the transmission probability T(e) and the Fermi functions /R,L(e) replaced 
by the effective transmission probability Tes{e) and the effective distribution 
function (89), respectively. 

Note that in general, ai^k 7^ o,N,k such that /R,eff 7^ /L,eff- This means that the 
driving can create an effective bias and thereby create a non-adiabatic pump 
current. Moreover, if all F„ are identical, the phase factors in (84) cancel each 
other and the effective Hamiltonian HeS equals the original static Hamiltonian. 

3.4-3 Linear-response limit 

For small driving amplitudes, it is often sufficient to treat the driving in the 
linear-response limit [146]. In doing so, we denote by g{t — t') the undriven 
limit of the Green function G{t,t') and by Hi{t) the time-dependent part of 
the Hamiltonian which is considered as a perturbation. Then, a formal solution 
of Eq. (33) is given by the Dyson equation 

/ + 00 
dt'g{t-t')ni{t')G{t',t-r), (90) 
-oo 

as can be shown by inserting (90) into (33). A self-consistent solution of this 
equation has been presented by Brandes [147]. Here, we restrict ourselves to 
the lowest order in the driving and, thus, can replace in the integral G{t',t — T) 
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by git' — t + t). Inserting moreover the Fourier representations 

ni{t)^ f^e-^-'n,{u;), (91) 

and Eq. (34), we obtain 

G{t,e)^g{e) + j ^e-'^' g{e + hw)n^{u) g{e). (93) 



For purely harmonic driving, 1-L\{t) — HiCOs{ftt), one finds for the Fourier 
coefficients (35) of the Green function the expressions 



G^^'\e) = lg{e±hn)H,g(e), 



(94) 
(95) 



while all Fourier components G^'^'^ with \k\ > 1, vanish to linear order. Con- 
sequently, the elastic transmission probability T^'^^e) is independent of the 
driving, i.e. it equals the result in the absence of external driving. The trans- 
mission probabilities under emission/absorption of a single photon are, how- 
ever, proportional to the intensity of the driving field, i.e. cx and read 



rL(e d= HQ) rR(e) {l\g{e ± HQ) Hi g{e)\N) 



(96) 



T^^\e) follows from the replacement (L, 1) ^ (R, A'"). 



3.5 Special cases 



In some special cases, the results of our Floquet approach reduce to simpler 
expressions. In particular, this is the case for zero driving amplitude, i.e. in 
the absence of driving, and for a driving that results from a time-dependent 
gate voltage and, thus, is homogeneous along the wire. 



3.5.1 Static conductor and adiabatic limit 

For consistency, the expressions (46) and (46) for the dc current and the zero- 
frequency noise, respectively, must coincide in the undriven limit with the cor- 
responding expressions of the time- independent scattering theory, Eqs. (17) 
and (21), respectively. This is indeed the case because the static situation is 
characterized by two relations: First, in the absence of spin-dependent interac- 
tions, we have time- reversal symmetry and therefore rLR(e) = TRL(e)- Second, 
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all sidebands with k vanish, i.e., Trl (^) — ^lr (^) — ^k,oT{e), where 

T(e) = rL(e)rR(e)|G'i^(e)|2 (97) 

and G{e) is the Green function in the absence of driving. Then the current as- 
sumes the known form (17). Moreover in a static situation, the matrix element 
of Eq. (38) reads [49] 

|rL(e)G'n(e)+i|' = l-r(e). (98) 

This relation allows one to eliminate the backscattering terms in the second 
line of Eq. (50) such that the zero-frequency noise becomes (21). Obviously 
if in a static situation both voltage and temperature are zero, not only the 
current (17) but also the noise (21) vanishes. In the presence of driving, this 
is no longer the case. This becomes particularly evident in the high-frequency 
limit studied in Section 3.4.2. 

It is known that in the adiabatic limit, i.e., for small driving frequencies, the 
numerical solution of the Floquet equation (55) becomes infcasiblc because a 
diverging number of sidebands has to be taken into account. In more math- 
ematical terms, Floquet theory has no proper limit as — > [148]. The 
practical consequence of this is that for low driving frequencies, it is favor- 
able to tackle the transport problem with a different strategy: If Ml is the 
smallest energy-scale of the Hamiltonian (1), one computes for the "frozen" 
Hamiltonian at each instance of time the current and the noise from the static 
expressions (46) and (50) being followed up by time-averaging. 

3.5.2 Spatially homogeneous driving 

In many experimental situations, the driving field acts as a time-dependent 
gate voltage, i.e., it merely shifts all on-site energies of the wire uniformly. 
Thus, the wire Hamiltonian is of the form 

n^ue(t)^no + f{t)J2\n){n\, (99) 

n 

where, without loss of generality, we restrict f{t) to possess zero time-average. 
A particular case of such a homogeneous driving is realized with a system that 
consists of only one level [118-120]. Then trivially, the time and the position 
dependence of the Floquet states factorize and, therefore, the dc current can 
be obtained within the formalism introduced by Tien and Gordon [117]. The 
corresponding noise properties have been addressed by Tucker and Feldman 
[115,116]. Here, we establish the relation between such a treatment and the 
present Floquet approach. 

Since the time-dependent part of the Hamiltonian is proportional to the unity 
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operator, the solution of the Floquet equation (55) is, besides a phase factor, 
given by the eigenfunctions \a) of the time- independent operator Tio — i^, 



ka(i))=e-^W|a), 
where {7io — iS)|Q!) = (e^ — ihja)\(^) and 



(100) 



(101) 



The quasienergies (e^ — i/i7a) coincide with the eigenvalues of the static eigen- 
value problem. Note that F{t) obeys the T-periodicity of the driving field 
since the time-average of f(t) vanishes. Thus, the phase factor in the Floquet 
states (100) can be written Fourier series, 



-iF(t) 



-ikUt 



(102) 



and, consequently we find \ua,k) = (ikloi) and the adjoint states accordingly. 
Then, the Green function (35) becomes 



G'^'He)=E«fc'+fc4'G'(e-A;'nQ), 



(103) 



where G{e) denotes the Green function in the absence of the driving field. 
Inserting (103) into (46) and employing the sum rule J2k' O't'^k'+k — Sk,o, 
yields 

^ = E / - khn)[Me) - /L(e)], (104) 

where ^(e) is the transmission probability in the absence of the driving. This 

expression allows the interpretation, that for homogeneous driving, the Flo- 
quet channels contribute independently to the current /. For the special case 
of a one-site conductor and a sinusoidal driving, this relation to the static 
situation has been discussed in Refs. [118, 119]. 

Addressing the noise properties, we obtain by inserting the Green function 
(103) into (50) the expression 



k ^ k' 



+ rLrR 

-|- same 



E 4'+feafe'G'ijv(e - k'hQ) [r^Gl^ie - k'Kl) - i 

k' 

terms with the replacement (L, 1) (R, A^)|. 



k{e)h{e + km) 



(105) 



While the term in the first line contains only the static transmission proba- 
bility at energies shifted by multiples of the photon energies, the contribution 
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in the second line cannot be brought into such a convenient form. The reason 
for this is that the sum over k' inhibits the application of relation (38). As 
a consequence, in clear contrast to the dc current, the zero-frequency noise 
cannot be interpreted in terms of independent Floquet channels. Only in the 
limit of large driving frequencies (cf. Section 3.4.2), the channels become ef- 
fectively independent and we end up with an expression that depends only 
on the transmission probability in the absence of the driving, and the Fourier 
coefficients a^. 

For large voltages where /l = and /r = 1, the sums over the Fourier 
coefficients in Eqs. (104) and (105) can be evaluated with the help of the 
sum rule J2k' o-t'^k'+k = ^fe,o- Then both the dc current and the zero-frequency 
noise become identical to their value in the absence of the driving. This means 
that for a transport voltage which is sufficiently large, a time-dependent gate 
voltage has no influence on the average current and the zero-frequency noise. 



4 Master equation approach 

An essential step in the derivation of the transmission within a weak-coupling 
approximation, Eq. (72), is the assumption that only terms with a — (3 and 
k = k' contribute significantly to (70). As discussed after Eq. (72), this requires 
that the separation of any pair of resonances is larger than their widths. This 
condition can be fulfilled only if the quasienergy spcctrTim docs not contain 
any degeneracies and if, in addition, the wire-lead coupling is very weak. Here, 
we refine the weak-coupling approximation scheme of Section 3.4.1 and derive 
a master equation approach which yields reliable results also in the presence 
of degeneracies and for intermediately strong wire-lead coupling [65,84]. 

4 ■ 1 Current formula 

We start again from the asymmetric expression (41) for the time-dependent 
current through the left contact. After averaging over the driving period, we 
obtain the dc current 



_ 9p f roo ri 

/^Ey dej^ dr <\tV,{e)U{e)l^^^^"^G\,{t,e)V^{T)Gu{t-T,e) 



for which we shall derive an approximation for small wire-lead coupling. 




(106) 
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Wc start with the second term which is hnear in the retarded Green function 
Gf'i{e). For small values of F, we obtain from (58) the approximation 

ImG^'\e) = 27r^ 6{e - e„ - kHQ) (107) 

a,k 

which allows one to perform the energy integration in Eq. (106). Then, we 
obtain the contribution 

-I E \ma,k)\' rL(e° + khQ) f{el + km - //l). (108) 

a,k 



The first term in Eq. (106) is quadratic in the Green function and, thus, 
requires a more elaborate treatment since otherwise, squares of 5-functions 
would emerge (cf. also the discussion in Section 3.4.1). For that purpose, it 
is advantageous to go one step back and to use instead of the current for- 
mula (106) the current operator (40) as a starting point. The time-average of 
the expectation value corresponding to the first term of Eq. (106) reads 



— / dr / di Re Ti^{T){c\{t)c^{t-T)) . (109) 
Til Jo Jo J 

Assuming that rL(e) is a slowly varying function in the relevant energy range, 
we can replace the time-evolution of Ci from the time t back to i — r by the 
Heisenberg operator ci{t — T,t) — C/o(i — r, t)ciUo{t — r, t) with Uq being the 
propagator (56) in the limit Fl/r — > 0. This means that Ci{t — r, t) represents 
the limit Fl/r of Ci(t — t). In order to include the coherent dynamics 
properly, it is convenient to introduce the "Floquet picture creation operators" 
Ca{t) which are defined by the time-dependent transformation [65,84] 

Ca(t)^Y.iMt)\n)cn. (110) 

n 

Using the inverse transformation c„ = Z]a(^l0a(^)) Ca(i), which follows from 
the completeness of the Floquet states at equal times, we obtain 

cn{t -T,t)^ Ee-''""* eH^°«+'='^W^(n|0«,,) c^{t). (Ill) 

a,k 

Inserting (111) with n = 1 into (109), we arrive at an expression that con- 
tains the time-dependent expectation values Pa/sit) = (c^(t) Ca{t))t with both 
operators taken at time t. The Papit) at asymptotic times, in turn, are deter- 
mined from a kinetic equation which we derive in the next subsection. Before 
doing so, however, we simplify Eq. (109) further by using of the fact, that at 
asymptotically long times, all Pafjif) become T-periodic functions and, thus, 
can be decomposed into a Fourier series Papif) = J2k^w{~^kflt)Pai3,k- This 
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brings Eq. (109) into the form 



- E 



dr Re 



r^ir) eH^°+^«^W'^(0;3,ifc+fe'|l)(l|0a,ik)Pa/3,fe'l ■ (112) 



By inserting for the lead response function r(r) its definition (28), we finaUy 
find for the time- averaged current through the wire the expression 



a,k ^(},k' 

-|(l|0„,fe)|V(ea + A;nQ-/.L) 



(113) 



Note that we have disregarded the principal value terms which correspond to 
an energy-renormalization due to the wire-lead coupling. 



^.2 Floquet- Markov master equation 



Having expressed the current in terms of the wire expectation values Pap{t), 
we now derive for them an equation of motion valid in the regime of weak 
to moderately strong wire-lead coupling. We thus consider the time-derivative 
Papit), which with the help of the zeroth-order Floquet equation (67), can be 
written as 



Pa^(i) = -^(6°-e°)P«Mi) + Tr 



p{t)cUt)co^{t) 



(114) 



For the evaluation of the second term on the right-hand side of the last equa- 
tion, we employ the standard master equation (23) presented in Section 2.5. 
Using twice the relation TrAlB^C] = Ti[A,B\C, which directly results from 
the cyclic invariance of the trace, we obtain 



Papit) 



1 /■°° / + ~ 
~ W Jo Ca{t), -^contacts], -ffcontacts(^ " T, t)] 



(115) 



For the further evaluation of Eq. (115), we write both i^contacts and i^coutacts ^ 
T,t) in terms of c„(t — T,t) for which we insert the approximation (111). 
After some algebra, we arrive at a closed differential equation for Papit). This 
determines the Fourier coefficients of the asymptotic solution, Pap,k, which 
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obey the inhomogeneous set of equations 



;° - e° - khn)P^p,k 



e=L,R k' 



+Te{e^^ + k'h^){(j)a,k'-k\ni){ni\(f)fs,k') /(e° + k'hQ - ne) 



- ^ei^a' + k''hfl){(j)a,k'+k''-k\rie){ne\(l)a',k'') Pa'p,k' 



- Y ^ei^p' + k"hQ){(j)fj'^k"\ne){ne\(j)fj^k+k"-k') Pai3',k' 

(3',k" 



(116) 



Here, we have assumed that the ideal leads always stay in thermal equilibrium 
and, thus, are described by the expectation values (7). Moreover, principal 
value terms stemming from an renormalization of the wire energies due to the 
couphng to the leads have again been neglected. 

The solution of the master equation (116) together with the current expres- 
sion (113) derived earlier, permits an efficient numerical calculation of the dc 
current through the molecular wire even for rather large systems or for energy- 
dependent couplings. Furthermore, as we shall exemplify below, this approach 
is still applicable in the presence of degeneracies in the quasienergy spectrum. 

4-3 Rotating-wave approximation 

The current formula (74) valid for very weak wire-lead coupling, which was 
derived in Section 3.4.1, can also be obtained from the master equation ap- 
proach within a rotating-wave approximation. Thereby, one assumes that the 
coherent oscillations of all Papit) are much faster than their decay. Then it is 
useful to factorize Pa/sit) into a rapidly oscillating part that takes the coherent 
dynamics into account and a slowly decaying prefactor. For the latter, one can 
derive a new master equation with oscillating coefficients. Under the assump- 
tion that the coherent and the dissipative time-scales are well separated, it 
is possible to replace the time-dependent coefficients by their time-average. 
The remaining master equation is generally of a simpler form than the orig- 
inal one. Because we work here already with a spectral decomposition of the 
master equation, we give the equivalent line of argumentation for the Fourier 
coefficients Pai3,k- 

It is clear from the Fourier representation of the master equation (116) that 



if 



- + khQ > Fl/r , 
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for ell a, (3,1, then the corresponding Pai3,k emerge to be small and, thus, may 
be neglected. Under the assumption that the wire-lead couplings are weak and 
that the Floquet spectrum has no degeneracies, the RWA condition (117) is 
well satisfied except for 

a^P, k^O, (118) 

i.e. when the prefactor of the l.h.s. of Eq. (116) vanishes exactly. This motivates 
the ansatz 

Pal3,k — Pa ^a,l3 h,0, (119) 

which has the physical interpretation that the stationary state consists of an 
incoherent population of the Floquet modes. The occupation probabilities Pa 
are found by inserting the ansatz (119) into the master equation (116) and 

E,.r,(60+™) ■ ^'"'^ 

Inserting this solution into expression (113) yields in the wide-band limit the 
current formula (74). 



4-4 Phonon damping 



In order to describe the electron transport under the influence of phonon 
damping, commonly a boson-like heat bath is coupled to each wire site, which 
renders the on-site energies fluctuating with quantum noise [21-23,26-32,66, 
149-154]. This can be considered as an extension of the spin-boson model to 
more than two sites and the presence of leads. For the master equation (23), 
one then has in the flrst line in addition the Hamiltonian of the phonon bath, 
while the electron-phonon coupling enters as a further dissipative contribution 
to the second hue. Note that this leaves the expression (113) for the current 
formally unchanged. 



4- 4-1 Hartree-Fock approximation 

When evaluating the master equation, however, it turns out that in addition 
to the terms containing the single-electron density matrix Paplt), two-electron 
expectation values of the form (cj cj^, C/3 c^)* appear. By iteration, one thus 
generates a hierarchy of equations up to A'"-electron expectation values. To 
obtain a description in terms of only the single-electron expectation values, 
we employ the Hartree-Fock decoupling scheme defined by the approximation 

(4 c\ Co) ^ (4 Co) {c\ Cfi) - (4 Cfi) {c\ Co) = PasPp^ " PfBsPa^- (121) 

Clearly, such a mean-field approximation only covers certain aspects of the 
full many-particle problem. Nevertheless, it offers a feasible and consistent de- 
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scription. As a most striking consequence, the Hartree-Fock decoupling (121) 
leaves the master equation non- linear [66]. 

4.4-2 Thermal equilibrium 

A potential problem of quantum master equations has been pointed out in 
Refs. [134,155], namely that they might not be consistent with the second 
law of thermodynamics — in particular, that they might not predict zero cur- 
rent even in the absence of both transport voltage and driving. This apparent 
lack of a proper equilibrium limit, however, is not inherent to master equa- 
tions of the form (23) themselves, but results from an inconsistent treatment 
at a later stage: It is crucial to employ in the second line of Eq. (23) the 
exact interaction picture operators of the uncoupled subsystems. Using any 
approximation indeed bears the danger of inconsistencies. Master equations 
whose equilibrium limit suffer from the mentioned problems, have, e.g. been 
derived in Rcf. [156] and applied to non-equilibrium situations with a finite 
transport voltage [93,94] and with time-dependent fields [150,157] where no 
contradiction occurs. 

Therefore, an important consistency check for quantum master equations is 
an equilibrium situation, where Hnn' is time-independent and where no ex- 
ternal bias is present (/x^ = n for all i). It can be demonstrated [66] that the 
final reduced master equation in the absence of both driving and voltage has 
the solution P^a = dadfa, with the population = f{Ea — /i), determined 
by the Fermi distribution and the energy Ea of the eigcnstates which 
represent the undriven limit of the Floquet states. Consequently, the current 
(113) vanishes in accordance with elementary principles of statistical physics. 



5 Resonant current-amplification 

A natural starting point for the experimental investigation of molecular con- 
duction under the influence of laser fields is the measurement of fingerprints of 
resonant excitations of electrons in the current-voltage characteristics. Treat- 
ing the driving as a perturbation, Keller et al. [146] have demonstrated that 
resonant electron excitations result in peaks of the current as a function of the 
driving frequency. Kohler et al. [67] included within a Floquet master equa- 
tion approach the driving exactly and later derived an analytical solution [68] 
which is in good agreement with an exact numerical solution. In a related 
work [158], Tikhonov et al. studied this problem within a so-called indepen- 
dent channel approximation [159] of a Floquet transport theory. As a central 
result, it has been found that, in particular for long wires, such excitations 
enhance the current significantly. In this section, we review the analytical 
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Fig. 3. Bridged molecular wire consisting of = 8 sites of which the first and the 
last site are coupled to leads with chemical potentials and fJ-R = hl + eV. 

treatment of Ref. [68] and compare to exact numerical results. 

As a working model we consider a so-called bridged molecular wire consisting 
of a donor and an acceptor site and N — 2 sites in between (cf. Fig. 3). Each of 
the A*" sites is coupled to its nearest neighbors by a hopping matrix element A. 
The dipole force (8) of the laser field renders each level oscillating in time with 
a position-dependent amplitude. The energies of the acceptor and the donor 
orbitals, |1) and |A^), are assumed to be close to the chemical potentials of the 
attached leads, Ei = ~ /xl ~ /^r- The bridge levels En, n = 2, . . . , N — 1, 
lie Eb ^ A, eV above the chemical potential. 



5.1 Static conductor 

Let us first discuss the static problem in the absence of the field, i.e. for A = 0. 
In the present case where the coupling between two neighboring sites is much 
weaker than the bridge energy, A <^ Eb, one finds two types of eigenstates: 
One group of states is located on the bridge. It consists of A^ — 2 levels with en- 
ergies in the range [Eb — 2A, Eb + 2A]. In the absence of the driving field, these 
bridge states mediate the super-exchange between the donor and the accep- 
tor. The two remaining states form a doublet whose states are approximately 
given by (|1) ± \N))/V2. Its splitting can be estimated in a perturbational 
approach [160] and is approximately given by 2A{A/Eb)'^^'^- Thus, the wire 
can be reduced to a two-level system with the effective tunnel matrix element 
Aba = Aexp{—K,{N — 2)), where k, = ln{EB/ A). If the chemical potentials 
of the leads are such that /iL > Eb and /iR < Ea, i.e., for a sufficiently large 
voltage, the current is dominated by the total transmission and for T ^ Aba 
can be evaluated to read 



— 



2e|Ap 

r 



e 



2k(7V-2) 



(122) 
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For the explicit calculation see, e.g., Ref. [91]. In particular, one finds an ex- 
ponentially decaying length dependence of the current [20.24, 161]. Moreover, 
in this limit, it is also possible to evaluate explicitly the zero-frequency noise 
to obtain the Fano factor F — S / e\I\ — 1. This value has a direct physical in- 
terpretation: Because the transmissions of electrons across a large barrier are 
rare and uncorrelated events, they obey Poisson statistics and, consequently, 
the mean number of transported electrons equals its variance. This translates 
to a Fano factor F = 1 [131]. 

5.2 Resonant excitations 

The magnitude of the current changes significantly when a driving field with 

a frequency Vt ^ Eb /h is switched on. Then the resonant bridge levels merge 
with the donor and the acceptor state to form a Floquet state. This opens a 
direct channel for the transport resulting in an enhancement of the electron 
current. 

In order to estimate the magnitude of the current through the resonantly 
driven wire, we disregard all bridge levels besides the one that is in reso- 
nance with the donor and the acceptor. Let us assume that this resonant 
bridge level [■0b) extends over the whole bridge such that it occupies the sites 
[2), . . . , |A^— 1) with equal probability — 2. Accordingly, the overlap be- 

tween the bridge level and the donor/acceptor becomes 

(lli^moleculelV's) = ^^'^^^^'^^ = -^^= = (^^^ | i^molecule | iV) (123) 

The resonance condition defines the energy of the bridge level as 
(V'sliJinoiecuielV's) = (recall that we have assumed Ejj — Ea — 0). 

We now apply an approximation scheme in the spirit of the one described in 
Ref. [91] and thereby derive a static effective Hamiltonian that describes the 
time- dependent system. We start out by a transformation with the unitary 
operator 

r A ^ 1 

S{t) = exp ] - i ^ \n){n\nt -ijT^Y. sin(ni) \. (124) 

Note that S{t) obeys the T-periodicity of the original driven wire Hamiltonian. 
As a consequence, the transformed wire Hamiltonian 

^molecule(i) = S\t)Hr.olecnle{t)S{t) - ihS\t)S{t) (125) 

is T-periodic as well. For hfl ^ A, we can separate time-scales and average 
-f^moiecuie(^) over the driving period. In the subspace spanned by |1), lips), and 
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|A^), the time-averaged wire Hamiltonian reads 



molecule,eff 



^dt 
r 



molecule 



^ 1 0^ 
1 1 



(126) 



with the effective tunnei matrix eiement 



and Ji the first-order Bessel function of the first kind. 



(127) 



The situation described by the Hamiltonian (126) is essentially the following: 
The central site \iPb) is coupled by matrix elements b to the donor and the 
acceptor site. Since the latter in turn couple to the external leads with a self 
energy r/2, their density of states is 



(128) 



Then, the tunneling of the electrons from and to the central site is essentially 
given by the golden rule rate 



(129) 



Like in the static case, we assume that the chemical potential of the left 
(right) lead lies above (below) the on-site energy of the acceptor (donor) and 
that therefore the donor is always occupied while the acceptor is always empty. 
Then, the rate of electrons tunnehng from the central site to the acceptor is 
given by the golden rule rate (129) times the occupation probability p of the 
state \iPb)- Accordingly, the rate of electrons from the donor to \iPb) is given 
by w times the probability 1 — p to find the central site empty. Consequently, 
the occupation of the resonant bridge level evolves according to 



p — w{l — p) — wp. 



(130) 



Equation (130) has the stationary solution p — 1/2 and, thus, for resonant 
excitations, the dc contribution of the time-dependent current is given by 



/res = ewp = e 



(131) 



Here, we have used for small arguments of the Bessel function the approxi- 
mation Ji{x) ~ X. The dc current (131) obeys an intriguing scaling behavior 
as a function of the wire length: Instead of the exponentially decaying length 
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Fig. 4. (a) Average current I as a function of the the driving frequency for 
various wire lengths A''. The scaled amplitude is ^ = O.IA; the apphed voltage 
I^R — I^L = 5A/e. The other parameters read F = O.IA and ksT = 0. (b) Average 
current for various driving ampUtudes A and couphng strengths F for a wire of 
length A^ = 8. (c) Fano factor F = S/el foi the wire length N = 8 and the 
wire-lead coupling F = O.IA. From Ref. [68]. 



dependence (122) that has been found for the static case, in the presence of 
resonant driving, a scaling I (x 1/N emerges. In particular for longer wires, 
this means that the external field enhances the conductance by several orders 
of magnitude. 
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5.3 Numerical results 



In order to corroborate the analytical estimates presented above, we treat the 
transport problem for the driven wire sketched in Fig. 3 numerically by solving 
the corresponding Floquet equation (55) and a subsequent evaluation of the 
expressions (46) and (50) for the dc current and the zero-frequency noise, re- 
spectively. For a wire with N = 5 sites, one finds peaks in the current when the 
driving frequency matches the energy difference between the donor/acceptor 
doublet and one of the — 2 = 3 bridge levels, cf. the sohd line in Fig. 4a. 
The applied voltage is always chosen so small that the bridge levels lie below 
the chemical potentials of the leads. In Figure 4a the scale of the abscissa is 
chosen proportional to {N — 1) such that it suggests a common envelope func- 
tion. Furthermore, we find from Fig. 4b that the dc current is proportional to 
A^/r provided that A is sufficiently small and F sufficiently large. Thus, the 
numerical results indicate that the height of the current peaks obeys 

W oc , (132) 

which is essentially in accordance with our analytical estimate (131). The main 
discrepancy comes from the fact that the overlap between the resonant level 
and the donor /acceptor differs from the estimate (123) by a numerical factor 
of the order one. Moreover, Fig. 4c demonstrates that at the resonances, the 
Fano factor assumes values considerably lower than one as expected for the 
transport through a resonant single level [52] . 



6 Ratchets and non-adiabatic pumps 

A widely studied phenomenon in driven transport is the so-termed ratchet 
effect: the conversion of ac forces without any net bias into directed mo- 
tion [69-74]. The investigation of this phenomenon has been triggered by the 
question whether an asymmetric device can act as a Maxwell demon, i.e., 
whether it is possible to ultimately convert noise into work. Feynman's fa- 
mous "ratchet and pawl" driven by random collisions with gas molecules, on 
first sight, indeed suggests that such a Maxwell demon exists. At thermal equi- 
librium, however, the whole nano-device obeys the same thermal fluctuations 
as the surrounding gas molecules. Therefore, consistent with the second law 
of thermodynamics, no directed motion occurs [162] and one has to conclude 
that the ratchet effect can be observed only in situations far from equilibrium. 

A basic model, which captures the essential physics of ratchets, is an asymmet- 
ric, periodic potential under the influence of an ac driving. In such a system. 
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even in the absence of any net bias, directed transport has been predicted for 
overdamped classical Brownian motion [69,72] and also for dissipative quan- 
tum Brownian motion in the incoherent regime [75, 76]. A related effect is 
found in the overdamped limit of dissipative tunneling in driven superlattices. 
There, the spatial symmetry is typically preserved and the directed trans- 
port is brought about by a driving field that includes higher harmonics of the 
driving frequency [163-165]. 

In the context of mesoscopic conduction, it has been found that the cyclic adia- 
batic change of the conductor parameters can induce a so-called pump current, 
where the charge pumped per cycle is determined by the area of parameter 
space enclosed during the cyclic evolution [42,44,166]. This relates the pump 
current to a Berry phase [43,77]. Beyond the adiabatic regime, pump effects 
have been investigated theoretically [46, 59, 78, 95, 167] and also been mea- 
sured in coupled quantum dots [38,93,168]. Since in the non-adiabatic regime, 
the main contribution to the pump current comes from electrons consider- 
ably below the Fermi surface, non-adiabatic electron pumping is essentially 
temperature independent [45]. 

The studies presented in this chapter were mainly motivated by two aspects: 

First, although infinitely extended "ideal" ratchets are convenient theoretical 
models, any experimental realization will have finite length, i.e., consist of a 
finite number of elementary units; cf. Fig. 5, below. Thus, finite size effects 
become relevant and it is intriguing to know the number of coupled wire units 
that are needed to mimic the behavior of a practically infinite system. Second, 
prior studies of quantum ratchets focussed on incoherent tunneling [75,76]. By 
contrast, the present setup allows one to investigate ratchet dynamics in the 
coherent quantum regime which has not been explored previously. 

The results of this section, have originally [65, 84] been computed for finite 
temperatures within the master equation approach of Section 4. In the limit 
of zero temperature, but otherwise equal parameters, the results from such a 
perturbative treatment agree essentially perfect with the corresponding exact 
solution obtained from Eq. (46). 



6. 1 Symmetry inhibition of ratchet currents 

It is known from the study of deterministically rocked periodic potentials [169] 
and of driven classical Brownian particles [170] that the symmetry of the 

equations of motion may rule out any non-zero average current at asymptotic 
times. Thus, before starting to compute ratchet currents, let us first analyze 
what kind of symmetries may prevent the effect. We consider situations, where 
the electron distributions in both leads are identical — in particular, situations 
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where both leads are in thermal equilibrium with a common chemical poten- 
tial, fhi^) = /r(c) = /(f) for all e. Then, no electromotive force acts and, 
consequently, in the absence of driving, all currents must vanish. An applied 
driving field, however, violates the equilibrium condition and can generate a 
finite dc current 



Obviously, the pump current vanishes if the condition Ti^^{e) = (^ ) 
fulfilled for all k and e. One might now ask whether this condition can be 
ensured by any symmetry relation. For the dipole driving considered here, 
the relevant symmetries are those studied in the Appendix A. 3, namely time- 
reversal symmetry, time-reversal parity, and generalized parity. In Section 3.3, 
we have already identified the symmetry-related channels which possess equal 
transmission probabilities. 

Looking at the relations (64), (65), and (66), it becomes clear that the gener- 
alized parity iSqp is the only symmetry that directly yields a vanishing pump 
current. This is so because it implies for the transmission probabilities the re- 
lation (66) and, thus, we find /pump = [65]. While time-reversal symmetry is 
without any consequence for the pump current, time-reversal parity has some 
rather subtle effect which follows from the fact that the transmission probabil- 
ities obey (65) and that in the weak-coupling limit of Section 3.4.1, in addition, 
relation (73) holds. Given these two relations, we obtain Tj^^lie) = T^^{e) and, 
thus, the dc current vanishes. Since the weak-coupling approximation is correct 
to lowest order in the coupling F, the consequence of time-reversal parity for 
quantum ratchets and Brownian motors is that we no longer find the generic 
behavior /pump oc F, but rather /pump oc F^. 

In the following, we consider two typical cases where the generalized parity 
symmetry is broken and, thus, a pump current emerges, namely (i) an asym- 
metric structure under the infiuence of a harmonic dipole force, the so-called 
rocking ratchet, and (ii) a spatially symmetric system for which generalized 
parity is broken dynamically by mixing with higher harmonics. 

6.2 Spatial symmetry-breaking: Coherent quantum ratchets 

A straightforward way to break generalized parity, is to use a conductor with 
an asymmetric level structure. Then, already a purely harmonic dipole driv- 
ing a{t) = Asin{Qt) in the Hamiltonian (8) is sufficient to generate a dc 
current. As a tight-binding model of such a structure, we have considered a 
wire consisting of a donor and an acceptor site and Ng asymmetric groups in 
the ratchet-hke configuration sketched in Fig. 5. In molecular structures, such 
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Fig. 5. Level structure of the wire ratchet with N = 8 sites, i.e., Ng = 2 asymmetric 
groups. The bridge levels are Eb above the donor and acceptor levels and are shifted 
by ±Es/2. 



an asymmetry can be achieved in many ways, and was explored as a source 
of molecular current rectification since the early work of Aviram and Rat- 
ner [3]. Later this effect has been found experimentally [7,99]. In general, an 
asymmetry can be created by attaching different chemical groups to the oppo- 
site sides of an otherwise symmetric molecular wire [7,99,171]. In our model, 
the inner wire states are arranged in A^^ groups of three, i.e. N — 2 = 3Ng. 
In each group, the first (last) level is lowered (raised) by an energy Es/2, 
forming an asymmetric saw-tooth like structure. The energies of the donor 
and the acceptor orbitals are assumed to be at the level of the chemical po- 
tentials of the attached leads and since no voltage is applied, we thus have 
El = En = /iL = f^R- The bridge levels En lie at Eb and Eb ± Es/2, as 
sketched in Fig. 5. 

Figure 6 shows the resulting stationary time-averaged current /. A quanti- 
tative analysis of a tight-binding model has demonstrated that the resulting 
currents lie in the range of 10^^ A and, thus, can be measured with today's 
techniques [84]. In the limit of very weak driving, we find / oc EgA'^ (Fig- 7). 
This behavior is expected from symmetry considerations: The asymptotic cur- 
rent must be independent of any initial phase of the driving field and therefore 
is an even function of the field amplitude A. This indicates that the ratchet 
effect can only be obtained from a treatment of the field beyond Kubo the- 
ory. For strong laser fields. Fig. 6 also shows that / is almost independent of 
the wire length. If the driving is intermediately strong, I depends in a short 
wire sensitively on the driving amplitude A and the number of asymmetric 
molecular groups Ng-. even the sign of the current may change with Ng, i.e. 
we find a current reversal as a function of the wire length. For long wires that 
comprise five or more wire units, the average current becomes again length- 
independent, as can be seen from Fig. 8. This identifies the observed current 
reversal as a finite size effect. As practical consequence, such relatively short 
wires can mimic the behavior of an (infinitely extended) quantum ratchet. 
Moreover, the fact that / converges to a finite value if the number of wire 
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Fig. 6. Time-averaged current through a molecular wire ratchet that consists of 
Ng bridge units as a function of the driving strength A. The bridge parameters are 
Eb = lOA, Es = A, the driving frequency is O = 3A/h, the couphng to the leads 
is chosen as Fl = Tr = 0.1 A/ h, and the temperature is ksT = 0.25 A. The arrows 
indicate the driving amplitudes used in Fig. 8. From Ref. [84]. 
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Fig. 7. Absolute value of the timc-avcraged current in a ratchet-like structure 
with Ng = 1 as a function of A'^Es demonstrating the proportionality to A^Es for 
small driving amplitudes. All other parameters are as in Fig. 6. At the dips on the 
right-hand side, the current / changes its sign. From Ref. [65] . 

units is enlarged, demonstrates that the dissipation caused by the coupling 
to the leads is sufficient to establish the ratchet effect in the limit of long 
wires. In this sense, no on-wire dissipation is required. Still, if the wire-lead 
model (1) is extended by electron- phonon coupling, the ratchet current might 
be enhanced [66]. 

Figure 9 depicts the average current vs. the driving frequency fi, exhibit- 
ing resonance peaks as a striking feature. Comparison with the quasienergy 
spectrum reveals that each peak corresponds to a non-linear resonance be- 
tween the donor/acceptor and a bridge orbital. While the broader peaks at 
hD, ~ Eb — lOA match the 1:1 resonance (i.e. the driving frequency equals 
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Fig. 8. Timc-averaged current as a function of the number of bridge units Ng for 
the driving amphtudes indicated in Fig. 6. The other parameters are as in Fig. 6. 
The connecting hnes serve as a guide to the eye. Prom Ref . [84] . 
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Fig. 9. Time-averaged current as a function of the driving frequency O for ^ = A 
and Ng = 1 (sohd hne). All other parameters are as in Fig. 6. The dotted line 
depicts the solution within the rotating- wave approximation (120). After Ref. [84]. 

the energy difference), one can identify the sharp peaks for HQ < 7A as 
multi-photon transitions. The appearance of these resonance peaks clearly 
demonstrates that the molecular bridge acts as a coherent quantum ratchet. 
Moreover, owing to the broken spatial symmetry of the wire, one expects an 
asymmetric current-voltage characteristic. This is indeed the case as depicted 
in Fig. 10. 

6.3 Temporal symmetry-breaking: Harmonic mixing 

The symmetry analysis in Section 6.1 explains that for a symmetric bridge 
without a ratchet-like structure as sketched in Fig. 3, the pump current (133) 
vanishes if the driving is a purely harmonic dipole force. This is so because 
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Fig. 10. Time-averaged current as a function of the applied static bias voltage V, 
which drops solely along the molecule. The driving amplitude is ^ = A, the driving 
frequency = 3A/h, and Ng = 1. All other parameters are as in Fig. 6. After 
Ref. [84]. 
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Fig. 11. Shape of the harmonic mixing field a{t) in Eq. (134) for Ai = 2A2 for 
different phase shifts 0. For cf) = 0, the field changes its sign for t ^ —t which 
amounts to the time-reversal parity »Stp- 

then the system is invariant under the generalized parity transformation Sqp 
and, thus, the transmission factors obey relation (66). Still, generalized parity 
can be broken in a dynamical way by adding a second harmonic to the driving 
field, i.e., a contribution with twice the fundamental frequency fl, such that 
it is of the form 

a{t) = Ai sm{Qt) + sin(2Qi + 0), (134) 

as sketched in Fig. 11. While now shifting the time t by a half period vr/l] 
changes the sign of the fundamental frequency contribution, the second har- 
monic is left unchanged. The generalized parity is therefore no longer present 
and we expect to find a non- vanishing average current. 
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Fig. 12. Average current response to the harmonic mixing signal with ampHtudes 
Ai = 2A2 = A, as a function of the couphng strength for different phase shifts 4>. 
The remaining parameters are O = 10A//1, = 5 A, k-^T = 0.25A, A*" = 10. The 
dotted hne is proportional to F; it represents a current which is proportional to F^. 
Prom Ref. [65]. 

The phase shift plays here a subtle role. For = (or equivalently any 
multiple of vr) the time-reversal parity iSxp is still present. Thus, according to 
the symmetry considerations in Section 6.1, the current vanishes within the 
weak-coupling approximation for the transmission probability, cf. Eq. (72). 
Since this approximation is only correct to linear order in F, the higher-order 
contributions typically remain finite and, consequently, for small coupling the 
pump current obeys I oc F^. Figure 12 confirms this prediction. Yet one ob- 
serves that already a small deviation from = is sufficient to restore the 
usual weak coupling behavior, namely a current which is proportional to the 
couphng strength F. This effect can be employed for the detection phase lags. 



Other features of the harmonic mixing current resemble the ones discussed 
above in the context of ratchet-like structures [65] . In particular, we again find 
for large driving amplitudes that the current becomes essentially independent 
of the wire length. Typically, the current reaches convergence for a length 



6.^ Phonon damping 



Including also the coupling of the wire electrons to a phonon heat bath, one 
can no longer employ the scattering formula (46) and for the computation of 

the dc current, one thus, has to resort to the master equation approach of 
Section 4. Here we only mention the main findings and refer the reader to the 
original work, Ref. [66]: The presence of phonon damping, generally increases 
the pump current up to one order of magnitude. This means that for quantum 
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ratchets, noise plays a rather constructive role. Moreover, phonon damping 
influences the dependence of the current on the phase lag by providing an 
additional shift towards a cos 4> behavior. 



7 Control setups 

A prominent example for the control of quantum dynamics is the so-called co- 
herent destruction of tunneling, i.e., the suppression of the tunneling dynamics 
of a particle in a double- well potential [85] and in a two- level system [85,87]. 
Recently, coherent destruction of tunneling has also been found for the dynam- 
ics of two interacting electrons in a double quantum dot [89,172]. A closely 
related phenomenon is the miniband collapse in ac-driven superlattices which 
yields a suppression of quantum diffusion [88,145,173]. In this chapter, we 
address the question whether a corresponding transport effect exists: If two 
leads are attached to the ends of a driven tunneling system, is the suppression 
of tunneling visible in conductance properties? Since timc-dcpcndcnt control 
schemes can be valuable in practice only if they operate at tolerable noise 
levels, the question is also whether the corresponding noise strength can be 
kept small or even be controlled. 

7. 1 Coherent destruction of tunneling 

In order to introduce the reader to the essentials of coherent destruction of 
tunneling in isolated quantum systems, we consider a single particle in a driven 
two-level system described by the Hamiltonian 

A A 

TiTLsit) = --2'^^ + 2"^^ cos(Qi). (135) 

If the energy of the quanta hQ of the driving field exceeds the energy scales 
of the wire, one can apply the high-frequency approximation scheme of Sec- 
tion 3.4.2 [87,91] and finds that the dynamics can be described approximately 
by the static effective Hamiltonian (84) which for the present case reads 

"^TLS.eflf = Y'^^' (^^^) 

with the tunnel matrix element renormalized according to 

A — > Aeff = Jo(A/hn)A. (137) 

Again, Jq denotes the zeroth order Bessel function of the first kind. If the 
ratio A/hfl equals a zero of the Bessel function Jq (i.e., for the values 2.405.., 
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Fig. 13. Level structure of the molecular wire with = 3 orbitals. The end sites 
are coupled to two leads with chemical potentials /xl and /xr = fJ-L — . 

5.520.., 8.654.., . . . ), the effective tunnel matrix vanishes and the tunnehng is 
brought to a standstill 

This reasoning is readily generalized to other tight-binding systems: If neigh- 
boring sites are coupled by a hopping matrix element A and the difference of 
their on-site energies oscillates with an amplitude A, one finds that the physics 
is determined by the renormalized matrix element (137), provided that h£l is 
the largest energy scale. 

1.2 Current and noise suppressions 

In order to investigate coherent destruction of tunneling in the context of 
transport, we consider the wire-lead setup sketched in Fig. 13 where the wire 
is described by the dipole Hamiltonian (8) with on-site energies En = 0. The 
wire is assumed to couple equally to both leads, Tl = Tr = T, and the 
numerical results are computed from the exact current formula (46). 

For a driven wire described by the Hamiltonian (8), it has been found [63,64,90] 
that the oscillating dipole force suppresses the transport if the ratio A/hQ 
is close to a zero of the Bessel function Jq- Moreover, in the vicinity of such 
suppressions, the shot noise characterized by the Fano factor (22) assumes two 
characteristic minima. These suppression effects are most pronounced in the 
high-frequency regime, i.e., if the energy quanta hfl of the driving exceed the 
energy scales of the wire. Thus, before going into a detailed discussion, we start 
with a qualitative description of the effect based on the static approximation 
for a high-frequency driving that has been derived in Section 3.4.2. 

Let us consider first the limit of a voltage which is so large that in Eq. (88), 
/R,ofr — /L,efr can be replaced by unity. Then, the average current is determined 
by the effective Hamiltonian 

N-l N 

Hes = -Acs E + |^+1)(^|) + E \n){n\, (138) 

n=l n=l 
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which has been derived by inserting the time-dependent part of the Hamil- 
tonian (8) into Eqs. (82) and (84). Then, obviously Tics is identical to the 
static part of the Hamiltonian (8) but with the tunnel matrix element renor- 
mahzed according to Eq. (137). Since the Bessel function Jq assumes values 
between zero and one, the amplitude of the driving field allows one to switch 
the absolute value of the effective hopping on the wire, Aes, between and A. 
Since the transmission probability of an undriven wire is proportional to |Ap, 
the effective transmission probability T'eff(e) acquires a factor Jq^A/HQ). This 
renormalization of the hopping then results in a current suppression [63,64,90]. 

For the discussion of the shot noise, we employ the Fano factor (22) as a 
measure. In the limit of large applied voltages, we have to distinguish two 
limits: (i) weak wire-lead coupling F <^ Acfr (i.e., weak with respect to the 
effective hopping) and (ii) strong wire- lead coupling F ^ A^s- In the first 
case, the tunnel contacts between the lead and the wire act as "bottlenecks" 
for the transport. In that sense they form barriers. Thus qualitatively, we face 
a double barrier situation and, consequently, expect the shot noise to exhibit 
a Fano factor F 1/2 [52]. In the second case, the links between the wire sites 
act as — 1 barriers. Correspondingly, the Fano factor assumes values F ^ 1 
ior N — 2 (single barrier) and Fail/2forA^ = 3 (double barrier) [174]. At 
the crossover between the two limits, the conductor is (almost) "barrier free" 
such that the Fano factor assumes its minimum. 

In order to be more quantitative, we evaluate the current and the zero- 
frequency noise in more detail thereby considering a finite voltage. This re- 
quires a closer look at the effective electron distribution (89); in particular, 
we have to quantify the concept of a "practically infinite" voltage. In a static 
situation, the voltage can be replaced by infinity, /R(e) = 1 = 1 — fh{^), if 
all eigenenergies of the wire lie well inside the range [/iLjAi^]. In contrast to 
the Fermi functions, the effective electron distribution (89) which is decisive 
here, decays over a broad range in multiple steps of size hQ. Since for our 
model, Teff(e) is peaked around e = 0, we replace here the effective electron 
distributions by their values for e = 0, 



for zero temperature. We have inserted the coefficients ai^k — Jk{A{N — 
l)/2Ml) and aN,k — J-kiA{N — l)/2h£t) which have been computed directly 
from their definition (86); Jk denotes the fcth order Bessel functions of the 
first kind. The current, the noise, and the Fano factor are given by the static 
expressions (17) and (21) with the transmission probability and the electron 
distribution replaced by the corresponding effective quantities, T^s and /es/, 




(139) 
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respectively. Thus, we obtain 



/=A/oo, (140) 



S=X^S^ + ^{l-Xy^, (141) 
F^XF^+^-^, (142) 

respectively, where the subscript oo denotes the corresponding quantities in 
the infinite voltage limit, 

/oo=^/dereff(e), (143) 

2 

Soo I dereff(e)[l - reff(e)], (144) 
and = Soo/eloo- The factor 

A = /R,eff(0)-/L,efr(0)= E ^li^^^^^^) (^^S) 

\k\<K{V) ^ ^"-^ ^ 

reflects the influence of a flnite voltage; K{y ) denotes the largest integer not 
exceeding e\V\/2Kl. Since Jk{x) ft; for |A;| > x and Y,k Jki^) ~ 1' 
A = 1 if K{V) > A{N — l)/2f\D,. This means that for small driving amplitudes 
A < eV/{N — 1), we can consider the voltage as practically inflnite. With an 
increasing driving strength, A decreases and, thus, the current becomes smaller 
by a factor A but still exhibits suppressions. By contrast, since < 1 for all 
situations considered here, we find from Eq. (142) that the Fano factor will 
increase with smaller A. 

Let us emphasize that unlike in the present case, the quenching of transmission 
observed in Refs. [175, 176] does not result from a renormalized inter-well 
tunnel matrix element, but rather originates from the appearance of the Bessel 
function Jq in the effective electron distribution (139). Therefore, at large 
voltages, this quenching will not give rise to current suppressions. 



7.3 Numerical results 



Figure 14a depicts the dc current and the zero-frequency noise for a wire with 
N — 3 sites and a relatively large apphed voltage, /xl — a*r = 50A. As a 
remarkable feature, we find that for certain values of the field amplitude A, 
the current drops to a value of some percent of the current in the absence 
of the field [63, 90] with a suppression factor which is fairly independent of 
the wire- lead coupling F [66]. The corresponding noise strength S exhibits 
similar suppressions and, in addition, has some small plateaus in the vicinity 
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Fig. 14. Time-averaged current / and zero-frequency noise S (a) as a function of 
the driving amplitude A for a wire with N = 3 sites with on-site energies En = 
and chemical potentials = —fih = 25A. The other parameters read J7 = 5A/h, 
r = 0.5A, and k^T = 0. Panel (b) displays the Fano F factor for these parameters 
(full line) and for smaller wire-lead coupling (dash-dotted line). From Ref. [63]. 

of the minima. The role of the plateaus is elucidated by the relative noise 
strength characterized by the Fano factor (22) which is shown in Fig. 14b. 
Interestingly enough, we find that the Fano factor as a function of the driving 
amplitude A possesses both a sharp maximum at each current suppression and 
two pronounced minima nearby. For a sufficiently large voltage, the Fano factor 
at the maximum assumes the value F ^ 1/2. Once the driving amplitude is 
of the order of the applied voltage, however, the Fano factor becomes much 
larger. The relative noise minima are distinct and provide a typical Fano factor 
of F 0.15. Reducing the coupling to the leads renders these phenomena 
even more pronoTinccd since then the suppressions occur in a smaller interval 
of the driving amplitude, cf. Fig. 14b. The overall behavior is robust in the 
sense that approximately the same values for the minima and the maximum 
are also found for larger wires, different driving frequencies, different coupling 
strengths, and slightly modified on-site energies, provided that A, F, En <^ Ml 
and that the appfied voltage is sufficiently large [64]. 

A comparison of these numerical results and the ones obtained in Section 7.2 
analytically within a high-frequency approximation shows an excellent agree- 
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Fig. 15. (a) Time-averaged current (solid line) and zero-frequency noise (dashed 
line) as a function of the driving amplitude for the driving frequency n = A/h and 
the transport voltage V = 48A/e. (b) Corresponding Fano factor for the same data 
(solid line) and for the driving frequencies = 1.5 A/ h (broken) and J7 = 3A/h 
(dash-dotted). All other parameters are as in Fig. 14. From Ref. [64]. 

merit: It quantitatively confirms both the parameter values for which current 
and noise suppressions occur and the corrections for in the large- amplitude 
regime A>eV [64,91]. 



For a much lower driving frequency of the order of the wire excitations, Q — 
A/H, the high-frequency approximation is no longer applicable. Nevertheless, 
the average current exhibits clear minima with a suppression factor of the order 
of 1/2; see Fig. 15a. Compared to the high-frequency case, these minima are 
shifted towards smaller driving amplitudes, i.e., they occur for ratios A/ Ml 
slightly below the zeros of the Bessel function Jq. At the minima of the current, 
the Fano factor still assumes a maximum with a value close to F 1/2 
(Fig. 15b). Although the sharp minima close to the current suppressions have 
vanished, in-between the maxima the Fano factor assumes remarkably low 
values oi F 0.2. 

So far, we have assumed that all on-site energies of the wire are identical. 
In an experimental setup, however, the applied transport voltage acts also a 
static dipole force which rearranges the charge distribution in the conductor 
and thereby causes an internal potential profile [112-114]. The self-consistent 
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Fig. 16. Time-averaged current (a) and Fano factor (b) as a function of the driving 
amplitude A for a wire with N = 3 sites in the presence of an internal bias. The 
on-site energies are Ei = b, E2 = 0, E3 = —b. All other parameters are as in Fig. 15. 
From Ref. [64]. 

treatment of such effects is, in particular in the time-dependent case, rather 
ambitious and beyond the scope of this work. Thus, here we only derive the 
consequences of a static bias without determining its shape from microscopic 
considerations. We assume a position-dependent static shift of the on-site en- 
ergies by an energy —bxn, i.e., for a wire with N — 3 sites. 

El = 6, E2^ 0, E3 = -b. (146) 

Figure 16a demonstrates that the behavior of the average current is fairly 
stable against the bias. In particular, we still find pronounced current sup- 
pressions. Note that since 6 ^ f2 a high-frequency approximation is still ap- 
plicable. As a main effect of the bias, we find reduced current maxima while 
the minima remain. By contrast, the minima of the Fano factor (Fig. 16b) 
become washed out: Once the bias becomes of the order of the wire-lead cou- 
pling, 6 ^ r, the structure in the Fano factor vanishes and we find F ^ 1/2 
for all driving amplitudes A < eV/{N — 1) [cf. the discussion discussion after 
Eq. (145)]. Interestingly, the value of the Fano factor at current suppressions 
is bias independent. 
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Fig. 17. Schematic top view of a setup where a molecule connected to three leads al- 
lows one to control the current flowing between the different leads (electro-chemical 
potentials fj,^, HCn P-Ci) ^ ^ function of the polarization angle i? of a linearly 
polarized laser field. 

7.4 Current routers 

So far, we have only considered driven transport through two-terminal de- 
vices. While the experimental realization of three and more molecular con- 
tacts is rather callenging, such systems can be described theoretically within 
the present formalism. As an example, wc consider a planar three-terminal 
geometry with = 4 sites as sketched in Fig. 17. Wc borrow from electrical 
engineering the designations E, Ci, and C2. Here, an external voltage is al- 
ways applied such that Ci and C2 have equal electro-chemical potential, i.e. 
/iCi — A*C2 7^ A*E- In a perfectly symmetric molecule, where all on-site energies 
are equal, reflection symmetry at the horizontal axis ensures that any current 
which enters at E is equally distributed among Ci_2, thus /d — Ic2 — ~-^e/2. 

The fact that this structure is essentially two-dimensional, brings about a new 
degree of freedom, namely the polarization of the laser field. We assume it to 
be linear with an polarization angle i9 as sketched in Fig. 17. The effective 
driving amplitudes of the orbitals that are attached to the leads acquire now 
a geometric factor which is only the same for both orbitals Ci and C2 when 
-3 — 0. For any other polarization angle, the mentioned symmetry is broken 
and the outgoing currents may be different from each other. The difference 
may be huge, as exemplified in Fig. 18. There, the current ratio varies from 
unity for ^9 = 0° up to the order of 100 for = 60°. Thus, adapting the 
polarization angle enables one to route the current towards the one or the 
other drain. 
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Fig. 18. Average currents (calculated within the master equation formalism) through 
contacts Ci (solid) and C2 (broken) as a function of the polarization angle ?? for 
the three-terminal device depicted in the Fig. 17. The chemical potentials are 
/7,E = ~/iCi = — = 50A; the on-site energies En = 0. The driving field is 
specified by the strength A = 25A and the angular frequency O, = IOA//I; the 
effective coupling is F = 0.1 A and the temperature ksT = 0.25A. From Ref. [90]. 

For a qualitative explanation of the mechanism behind this effect, it is instruc- 
tive to look at the time-averages of the overlaps \{n\(f)a{t))\'^ = J2k K^l^a.fe)^ 
of the Floquet states with the terminal sites n — E, Ci, C2, which determine 
the effective tunneling rates (75) and (76) in the weak wire-lead coupling limit. 
Figure 19 shows these overlaps for three different polarization angles i). Let us 
consider, for instance, the current across contact Ci. It is plausible that only 
Floquet modes which have substantial overlap with both the site Ci and also 
the site E contribute the current through these terminals. For a polarization 
angle ^9 = —60°, we can infer from Fig. 19 that the Floquet states with indices 
a = 1, 3 and 4 fulfill this condition and, consequently, a current flows from 
lead E into lead Ci. By contrast, for = 0° and = 60° such current carrying 
states do not exist and the respective current vanishes. 



7.5 Phonon damping 



A further question to be addressed is the robustness of the current suppressions 
against dissipation. In the corresponding tunneling problem, the driving alters 
both the coherent and the dissipative time scale by the same factor [178]. 
Thus, one might speculate that a vibrational couphng leaves the effect of the 
driving on the current qualitatively unchanged. This, however, is not the case: 
With increasing dissipation strength, the characteristic current suppressions 
become washed out until they finally disappear when the damping strength 
becomes of the order of the tunnel coupling A [66] . This detracting influence 
underlines the importance of quantum coherence for the observation of those 
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Fig. 19. Time-average ((0a(t)|n)(n|^a(i))) of the overlaps \{n\4>a(t))^ of the sites 
n = E,Ci,C2, and T (central site) to a Floquet state \(j)ci{t)) for three different 
polarization angles 'd. All parameters are as in Fig. 18. Adapted from Ref. [177] 



current suppressions. Moreover, for the model employed in Ref. [66], we do 
not find the analogue of the effect of a stabilization of coherent destruction 
of tunneling within a certain temperature range [179-181] or, likewise, with 
increasing external noise [182], as it has been reported for driven, dissipative 
symmetric bistable systems. 



8 Conclusion and outlook 



In the present survey we have reviewed the role of external driving for various 
transport quantities in nanosystems. In particular, we have focussed on the 
possibilities to selectively control, manipulate and optimize transport through 
such systems. In this context, we have studied various aspects of the electron 
transport through time-dependent tight-binding systems. For the theoretical 
description, two formalisms have been employed which both take advantage 
of the Floquet theorem: A Floquet scattering approach provides an exact 
solution of the time-averaged electrical current beyond the linear response 
limit and, moreover, yields an expression for the corresponding noise power. 
Interestingly, unlike in the undriven case, the noise depends also on the phases 
of the transmission amplitudes. As a drawback, this scattering approach is 
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limited to the case of purely coherent transport in the absence of electron- 
electron interactions. As soon as other degrees of freedom like, e.g. a phonon 
bath, come into play, it is advantageous to resort to other formalisms like a 
Floquet master equation approach which, however, is limited to a weak wire- 
lead coupling. 

We have investigated several driven transport phenomena such as resonant 
current amplification (Section 5), electron pumping (Section 6), and coherent 
current control (Section 7). Of foremost interest in view of ongoing experi- 
ments is the enhancement of molecular conduction by resonant excitations. 
We have derived an analytical expression for the current enhancement factor 
and, moreover, have found that the relative current noise is reduced approxi- 
mately by a factor of one half. 

Both molecular wires and quantum dot arrays can act as coherent quantum 

ratchets and thereby operate in a regime where the quantum ratchet dynamics 
has not been studied previously. Of particular practical relevance is the fact 
that already relatively short wires or arrays behave like infinite systems. For 
the investigation of such driven nano-devices, symmetries play a crucial role: 
The driven nano-system may exhibit a dynamical symmetry which includes a 
time transformation. Breaking this dynamical symmetry, for instance by using 
a non-harmonic driving force, can be exploited for the generation of a pump 
current. Moreover, the symmetry analysis revealed that a ratchet or pump can 
only be observed in the absence of the so-called generalized parity. 

Coherent destruction of tunneling has a corresponding transport effect which 
exhibits an even richer variety of phenomena. For driving parameters, where 
the tunneling in isolated unbiased systems is suppressed, the dc current drops 
to a small residual value. This efi'ect is found to be stable against a static bias. 
Moreover, the investigation of the corresponding noise level characterized by 
the Fano factor, has revealed that the current suppressions are accompanied 
by a noise maximum and two remarkably low minima. This allows one to se- 
lectively control both the current and its noise by ac fields. Of crucial interest 
for potential applications are the noise properties of non-adiabatic pumps. For 
resonant excitations, these can be treated analytically within an approxima- 
tion scheme in the spirit of the one apphed in Ref. [91]. 

An experimental realization of the phenomena discussed in this paper is obvi- 
ously not a simple problem. The requirement for asymmetric molecular struc- 
tures is easily realized as discussed above, however difficulties associated with 
the many possible effects of junction illumination have to be surmounted [183]. 
Firstly, there is the issue of bringing the light into the junction. This is a diffi- 
cult problem in a break-j unction setup but possible in an scanning probe mi- 
croscope configuration. Secondly, in addition to the modulation of electronic 
states on the molecular bridge as discussed in this work, other processes involv- 
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ing the excitation of the metal surface may also affect electron transport. A 
complete theory of illuminated molecular junctions should consider such pos- 
sible effects. Moreover, the junction response to an oscillating electromagnetic 
field may involve displacement currents associated with the junction capacity. 
Finally, junction heating may constitute a severe problem when strong electro- 
magnetic fields are applied. On the other hand, the light-induced rectification 
discussed in this paper is generic in the sense that it does not require a par- 
ticular molecular electronic structure as long as an asymmetry is present. The 
prediction that with proper illumination, one might induce a unidirectional 
current without any applied voltage, offers the possibility to observe a pump 
effect without the background of a voltage-induced direct current component. 

An alternative experimental realization of the presented results is possible 
in semiconductor hetero-structures, where, instead of a molecule, coherently 
coupled quantum dots [37] form the central system. A suitable radiation 
source that matches the frequency scales in this case must operate in the 
microwave spectral range. Compared to molecular wires, these systems by 
are well-established. This is evident from the fact that in microwave- driven 
coupled quantum dots, electron pumping has already been observed [38]. 

The authors share the behef that this survey on driven quantum transport on 
the nanoscale provides the reader with a good starting point for future own 
research: Many other intriguing phenomena await becoming unraveled. 
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A A primer to Floquet theory 



In this review, we deal with time-periodically driven quantum systems whose 
dynamics is governed by the Schrodinger-like equation of motion 

ih^^m))-{H{t)-i^)\m) (A.l) 

with the T-periodic Hamiltonian H{t) = H{t + T). The hermitian self-energy 
term E results from an elimination of environmental degrees of freedom and 
renders the time-evolution non-unitary 

The explicit time-dependence in the Hamiltonian rules out the standard sepa- 
ration ansatz \ip{t)) = exp{—iEt/h)\ip), where E is the (complex) eigenenergy 
of a state for the solution of Eq. (A.l). Yet, the time-periodicity of the 
Hamiltonian allows one to apply Floquet theory, a powerful tool, which we 
briefly review in this appendix. 

A.l Floquet theorem for non-unitary time-evolution 

Floquet theory is based on the Floquet theorem which states that for a time- 
periodic Hamiltonian, H{t) = H{t + T), there exists a complete set {\'ipa{t))} 
of solutions of Eq. (A.l) which is of the form 

\Mt)) = e-('^«/'^+^")* \u^{t)) , \ua{t)) - Mt + T)) . (A.2) 

The time-periodic functions \ua{t)) are called Floquet modes or Floquet states 
and the quantities are referred to as quasienergies with corresponding 
width 7q. By inserting the ansatz (A.2) into Eq. (A.l), one easily verifies 
that the Floquet states fulfill the eigenvalue equation 

[H{t) -iE-ih ^) Kit)) = (e« - ih'ya) \ua{t)) . (A.3) 

Different methods can be used to prove the Floquet theorem. Here, we present 
a constructive argument. Upon diagonalization of the one-period propagator 
U{T,0), where U{t,t') is the time-evolution operator corresponding to the 
dynamical equation (A.l), we obtain 

C/(r, 0) lu^m = e-(-«/^+^«)^ |^„(0)) . (A.4) 

Here, we have written the complex eigenvalue as exponential for some Cq, and 
ja- Next, we use the eigenstates \ua{0)) as initial states for the time-evolution 
according to Eq. (A.l), yielding the solutions li'ait)) = U{t,0)\ua{0)) of 
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Eq. (A.l). This allows us to define the Floquet modes \ua{t)) — exp[{iea/h + 
la)'t]\'4^ot{'t)) , which are indeed T-periodic functions: 

\u^{t + T)) = e('^«/'^+^«)(*+^) U{t + r, 0)|ii«(0)) 

= e('^«/'^+^«)(*+^) U{t, 0) C/(r, 0) |iic«(0)) (A.5) 
= eP^«/'^+^«)* iVa(t)) = l^a(i)) ■ 

In the second line, we have used that owing to the time-periodicity of the 
Hamiltonian, the relation U{t + T,T) = U{t, 0) holds true for arbitrary time t. 
Finally, the completeness of the set of solutions {|'0a(^))} follows, if we assume 
the completeness of the eigenstates of U (T, 0) . 

Since the one-period propagator U (T, 0) is in general non-unitary, its eigen- 
states \ua{0)) are not mutually orthogonal. We therefore also have to consider 
the left eigenstates of C/(T, 0), i.e., the solutions of the adjoint equation 

I^H{t) + iE^-m^y+{t)) = {e^ + ih^^)K{t)) . (A.6) 

Here, we have used the fact that the eigenvalues of the adjoint equation are 
the complex conjugates of the eigenvalues of the original eigenvalue equation 
(A. 3). This follows from the secular equations corresponding to the eigen- 
value problems (A. 3) and (A.6) by using the relation det O = det O^, which 
holds for an arbitrary operator O. Assuming completeness of the eigenstates 
of C/(T, 0), the Floquet modes and its adjoint modes may be chosen to form 
a bi-orthonormal basis at equal times t, 

{uM\Mt)) = ^ap and Y.\<it))Mt)\ = l- (A.7) 

a 

The time-evolution operator U{t, t') can be expressed explicitly in terms of 
the Floquet modes and quasi-energies to read 

U{t,t') = J2e-'^'"/^+'r.){t-t') \^^^t)){u+{t')\ . (A.8) 

a 

This relation is readily checked by noting that due to Eq. (A. 2) the right-hand 
side solves the differential equation (A.l). The initial condition U{t,t) = 1 is 
ensured by the completeness (A.7) of the Floquet modes. 

It is worthwhile to remark that the conceptual importance of Floquet the- 
ory lies in the fact that it allows one to separate the long-time dynamics, 
governed by the eigenvalues — ih'ja, from the dynamics within one driving 
period, determined by Floquet modes \ua{t)) [cf. Eq. (A. 2)]. Note also that 
the quasienergies and the Floquet states in Eq. (A. 2) are not defined uniquely. 
In fact, the replacement 

^ ea + Khn , \u„{t)) ^ e^'^"^* \Ua{t)) , (A.9) 
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where ft — 27r/T is the angular frequency of the driving and {ka} is an 
arbitrary sequence of integers, yields a new set of quasienergies and Floquet 
states corresponding to the same solutions Hipait))} of Eq. (A.l). In other 
words, the quasienergies and Floquet modes come in classes, out of which one 
is allowed to select a single representative, usually with quasienergy in a single 
"Brillouin zone" E — hQ/2 < ea < E + hQ/2, where E is an arbitrary but 
fixed energy. 



A. 2 Extended Hilbert space formalism 



According to the basic postulates of quantum mechanics, the state of a system 
is described by a vector l'^) in a Hilbert space R with the inner product 
Without loss of the generality, we assume that there exists a countable and 
complete set {\n)} of orthonormal states, i.e., 

{n\n') = 5nn' , Yl = 1 ■ 

n 



The Hilbert space T of all T-periodic, complex-valued functions possesses the 
inner product 

1 



{u,v) = - / dtu*it)vit) (A.ll) 
i Jo 

and the functions exp{ikQt) with k = 0, ±1, ±2, . . . form the corresponding 
complete and orthonormal set. The decomposition of an arbitrary T-periodic, 
complex-valued function into this basis yields the standard Fourier series. 

As first noted by Sambe [138], the time-periodicity of the Floquet modes 
suggests their description in the composite Hilbert space M (g) T. Its elements, 
for which we adopt the notation \u)) [138], are the T-periodic state vectors 
\u{t)) = \u{t+T)). Introducing the inner product in this space in the canonical 
way via 

{{u'\u)) = ^l^dt{u'(tMt)) , (A.12) 
an orthogonal basis of R T is given by the set of states defined by 

\u^^{t)) = e'*^"* \n) . (A.13) 

The arbitrary integer k is sometimes called the sideband index. The decompo- 
sition of a state \u{t)) into this basis is equivalent to the Fourier representation 



k 

\ua,k) = ^ e''''' \u^it)) 
1 Jo 



(A.14) 
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It is important, however, to always keep in mind that the states \ua,k) are not 
orthogonal, because the Floquet modes are only mutually orthogonal at equal 
times [cf. Eq. (A. 7)]. 

By the introduction of a Hilbert space structure for the time dependence, we 
have formally traced back the computation of Floquet states to the compu- 
tation of eigenstates of a time-independent Hamiltonian with an additional 
degree of freedom. In particular, in the composite Hilbert space the Floquet 
equation (A. 3) maps to the time- independent eigenvalue problem 

{n{t)-iT.)\u^))^e^\u^)) (A.15) 

with the so-called Floquet Hamiltonian 

n{t) = H{t) - ihj^ . (A.16) 



A wealth of methods for the solution of this eigenvalue problem can be found 
in the hterature [139, 184]. One such method is given by the direct numerical 
diagonalization of the operator on left-hand side of Eq. (A.15). For a har- 
monic driving, the eigenvalue problem (A.15) is band-diagonal and selected 
eigenvalues and eigenvectors can be computed by a matrix-continued fraction 
scheme [184,185]. 

In cases where many Fourier coefficients (in the present context frequently 
called "sidebands") must be taken into account for the decomposition (A. 14), 
direct diagonalization is often not very efficient and one has to apply more 
elaborated schemes. For example, in the case of a large driving amplitude, one 
can treat the static part of the Hamiltonian as a perturbation [87,138,145]. 
The Floquet states of the oscillating part of the Hamiltonian then form an 
adapted basis set for a subsequently more efficient numerical diagonalization. 

A completely different strategy to obtain the Floquet states is to propagate 
the Schrodinger equation for a complete set of initial conditions over one 
driving period to yield the one-period propagator. Its eigenvalues represent 
the Floquet states at time t = 0, i.e., \ua{Q)). Fourier transformation of their 
time-evolution results in the desired sidebands. Yet another, very efficient 
propagation scheme is the so-called (i, i')-formalism [186]. 



A. 3 Parity of a system under dipole driving 



Although we focus in this work on tight-binding systems, it is more convenient 
to study symmetries as a function of a continuous position and to regard the 
discrete models as a limiting case. Moreover, we consider in this section the 
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Hamiltonian of the entire system including the leads. Consequently, we do not 
have not include any self-energy contribution. 

A static Hamiltonian Hq{x) is called invariant under the parity transformation 
V : X ^ —X if it is an even function of x. Then, its eigenfunctions can be 
divided into two classes: even and odd ones, according to the sign in ipaix) = 
=t<^a(~2;). Adding a periodically time-dependent dipole force xa{t) to such a 
Hamiltonian evidently breaks parity symmetry since V changes the sign of the 
interaction with the radiation. In a Floquet description, however, we deal with 
states that are functions of both position and time — we work in the extended 
space ]R ® T. Instead of the stationary Schrodinger equation, we address the 
eigenvalue problem 

n{x, t) (j){x, t) = e (l){x, t) (A. 17) 

with the Floquet Hamiltonian for zero self-energy given by 

n{t) = Ho(x) + xa(t) - ih^, (A.18) 

ot 

where we assume a symmetric static part, Hq{x) = Hq{—x). Our aim is now 
to generalize the notion of parity to the extended space M (g) T such that the 
overall transformation leaves the Floquet equation (A. 17) invariant. This can 
be achieved if the shape of the driving a{t) is such that an additional time 
transformation "repairs" the acquired minus sign. We consider two types of 
transformation: generalized parity and time-reversal parity. Both occur for 
purely harmonic driving, a{t) = cos{Qt). In the following we derive their 
consequences for the Fourier coefficients 

M^)-^ / di e'^"^ i) (A.19) 
I Jo 

of a Floquet states (f){x,t). 



A. 3.1 Time-reversal symmetry 

Before discussing parity symmetry, let us comment on time-reversal symmetry 
which is not relevant for the spectral properties but still has some computa- 
tional importance. It is known that the energy eigenfunctions of an non-driven 
Hamiltonian, which obeys time- reversal symmetry, can be chosen real [143]. 
Time-reversal symmetry is typically broken by a magnetic field (recall that 
a magnetic field is described by an axial vector and, thus, changes its sign 
under time-reversion) or by an explicit time-dependence of the Hamiltonian. 
However, for a cosine driving, time-reversal symmetry 

St: -t, (A.20) 

is retained and the Floquet Hamiltonian (A.18) obeys 'H{t) = [H{—t)]*. With 
the same line of reasoning as in the case of time-reversal symmetry, but with 



69 



the additional replacement x — > — x, we obtain that one can choose the Floquet 
states with (j){x,t) — (j)*{x, —t). Then, the Fourier coefficients (A. 19) are real 

= <Pl{x), (A.21) 

which helps to reduce numerical effort. 

A. 3. 2 Time-reversal parity 

A further symmetry is found if a is an odd function of time, a(t) = —a{—t), e.g. 
for a{t) = sm[Qt). Then, time inversion transforms the Floquet Hamiltonian 
(A. 18) into its complex conjugate such that the corresponding symmetry is 
given by the anti-linear transformation 

cStp : {(p,x,t) ^ {(P*,-x,-t). (A.22) 

This transformation represents a generalization of the parity V; we will refer 
to it as time-reversal parity since in the literature the term generalized parity 
is mostly used in the context of the transformation (A. 24). 

Again we are interested in the Fourier decomposition (A. 19) and obtain 

0,(x)=0*(-x). (A.23) 

The time-reversal discussed here can be generalized by an additional time- 
shift to read t ^ — t. Then, we find by the same line of argumentation that 

(f)k{,x) and (pki—x) differ at most by a phase factor. However, for convenience 
one may choose already from the start the origin of the time axis such that 
to = 0. 

A. 3. 3 Generalized parity 

It has been noted [85,86, 144] that a Floquet Hamiltonian of the form (A. 18) 
with a{t) — sin(Qi) may possess degenerate quasienergies due to its symmetry 
under the so-called generalized parity transformation 

Sgp ■■ {x, t) {-X, t + Ti/Q), (A.24) 

which consists of spatial parity plus a time shift by half a driving period. 
This symmetry is present in the Floquet Hamiltonian (A. 18), if the driving 
field obeys a(t) = —a{t + vr/f2), since then Sqp leaves the Floquet equation 
invariant. Owing to S^p = 1, we find that the corresponding Floquet states 
are either even or odd, Sgp4>{x, t) = (j){—x,t-\-7i/fl) = ±(j){x, t). Consequently, 
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the Fourier coefficients (A. 19) obey the relation 



(A.25) 
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